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Large Dislocation Loops in Antimony Telluride} 


By P. DELavicNETTE and S. AmELINCKx 
Centre d’Etude de I’Energie Nucleaire, MOL, Belgium 


[Received November 17, 1960] 


ABSTRACT 


When slowly grown from the melt, antimony telluride contains large 
hexagonal loops. The observations suggest that these are prismatic loops 
resulting from deviations of the stoichiometric composition. In particular 
it is proposed that tellurium vacancies precipitate in prismatic loops, from 
which the ‘ geometrical’ stacking fault is eliminated, by glide of a partial. 
An explanation is given for the occurrence of pairs of concentric loops 
elongated in direction differing by 60°. The equilibrium shape as calculated 
on the basis of the model agrees reasonably well with the observed shape. 


§ 1. INTRODUCTION 


In the course of a study of dislocation in Sb, Te, an interesting precipitation 
phenomenon was found. A brief description was already presented at 
the Delft Meeting on Electron Microscopy. The phenomenon is as yet 
not completely understood but we believe it worth while to present the 
experimental facts and suggest a possible explanation. These observations 
may possibly provide some insight in the question how stoichiometric 
deviations are accommodated in a crystal. 


§ 2, EXPERIMENTAL 


Crystals of antimony telluride were prepared by melting together the 
elements in the stoichiometric proportion in evacuated and sealed quartz 
capsules. The melt was cooled very slowly in order to obtain large grains. 
Specimens for electron microscopic observation were prepared in the 
way described earlier, i.e. by cleavage (Delavignette and Amelinckx 1960). 
It was found that the specimens contained invariably large elongated 
hexagonal loops of the kind shown in fig. 1(c){. The lines exhibit diffraction 
contrast. We will further present evidence that they are prismatic disloca- 
tion loops. 

In more rapidly cooled specimens smaller rounded loops were found 
(fig. 1 (a,b)); in quenched specimens no loops at all were observed suggesting 
that they are not due to thermally generated frozen-in point defects. In 
other words the loops presented here are not quench loops. Qualitative 
experiments seem to indicate that the number and total area covered by 
loops depend on the composition of the melt, from which the crystal grew. 


+ Communicated by Dr. J. Goens. 
{ Figures 1 and 2 are shown as plates. 
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In specimens which were very slowly cooled very often concentric 
pairs of large loops are found. The two loops are then systematically 
elongated in directions which differ by 60°, as shown in fig. 2. The sides of 
the hexagons are in all cases parallel to the close-packed directions of the 
lattice. Other dislocations interact visibly with the hexagons, especially 
with the small ones (fig. 1(6)). The loops are sessile; whereas other dis- 
locations frequently move, no such movement was ever noticed for the 
loops. Bismuth telluride specimens treated in a similar way did not 
show analogous features. 


§ 3. NATURE OF THE LOOPS 


From most of the photographs it is evident that the loops show an 
inhomogeneous contrast; certain segments give practically no contrast. 
If the set of planes that gives rise to the contrast is perpendicular to the 
short segment, no contrast is observed for that segment, the others show 
some. This is represented schematically in fig. 3 (a, 6). If the set of 
diffracting planes is perpendicular to one of the larger segments all segments 
give contrast. These observations are in accord with the theory of Howie 
and Whelan (1961) on the contrast from prismatic dislocation loops. The 
observations further suggest that the Burgers vector should be inclined 
with respect to the plane of the loop, in the direction perpendicular to the 
short segment. 

No stacking fault contrast is found in the loops. 


Fig. 3 


(a) (b) 


Schematic representation of the contrast effects observed at loops. The bar-dot 
line represents the trace of the reflecting planes. 


§ 4, Discussion 


4.1. The Crystal Structure of SbyTes 


The crystal structure of Sb,Te, is isomorphous with that of Bi,Te,. 


It consists of a cubic stacking of close-packed layers of antimony and 
tellurium, 
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If greek letters represent antimony and latin letters tellurium. the 
stacking sequence in one unit cell can be symbolized as: 


aBcub, cabya, byaBe, ... 


where a, (or «), 6 (or 8) and ¢ (or y) stand for the three possible positions. 
It is clear that tellurium occurs in two different positions, we call Te, 
the atom sites limiting the five-layered lamellae, Te, the sites in the centre 
of these lamellae. The bonds between the two Te, layers are the weakest 
and the glide plane is therefore very probably located there. 

Antimony telluride prepared from the melt seems to be invariably 
non-stoichiometric, in particular containing an excess of antimony 
(Offergeld and Van Caekenberghe 1959). It is further generally thought 
that the disorder in this compound is of the anti structure type, i.e. antimony 
occupying tellurium sites and vice versa (Kroger and Vink 1956). 


4.2. Generation Mechanism of the Loops 

Although the origin of the loops is not completely understood it is fairly 
clear that the reason for their occurrence has to be sought in the deviations 
from stoichiometry. This point of view is based on the fact that similar 
features do not occur in bismuth telluride (which is nearer to stoichiometry 
when grown from the melt) and also because the number and size of the 
loops were found to be sensitive to deviations from stoichiometry. It 
should be mentioned that these deviations were in all cases smaller than 
those needed to give rise to a second phase. 

It is now clear that the crystal should contain a number of point defects 
which make up for the non-stoichiometry. In view of the type of dis- 
order, one would expect at high temperatures more antimony atoms on 
tellurium sites than vice versa. The vacancies cannot be given a name as 
a consequence of the disorder, since they may occupy both types of sites. 
By cooling slowly the atoms will tend to occupy their proper sites and 
this will result in a preferential occupation of the Te layers by vacancies. 
Since bonds are weakest and diffusion easiest between two Te,-layers, 
vacancy loops will preferentially nucleate in one of these. There is a 
further reason why this type of layer is preferred, but this will only become 
apparent a posterio. 

The type of loops that would result is shown in fig. 4(a). It contains 
a ‘chemical’ and a ‘geometrical’ fault. Since the basal dislocations in 
these crystals are not dissociated the stacking-fault energy is relatively 
large and the geometrical stacking fault will therefore be eliminated 
out of the loop by the nucleation of a partial, as in aluminium (Hirsch 
et al. 1958). Since the plane of the loop is a glide-plane this should 
be not too difficult; it would be difficult if the loop were to nucleate 
in another layer than a Te;-layer. The final result is then as shown in 
fig. 4(b). The Burgers vector is now inclined with respect to the c-plane 
in agreement with the observed contrast effect. The two pure edge 
segments presumably climb somewhat more rapidly than the four segments 


2Q2 
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of mixed character. In this way a natural explanation for the elongated 
shape is found. We show further that the equilibrium shape is such that 
the edge segments should be shorter. 


Fig. 4 
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(¢) 


Arrangement of layers in the prismatic loops. (a) Loop without shift; (b) 
loop after elimination of geometrical stacking fault; (c) two concentric 
loops. 


The loops do not contain a geometrical stacking fault, i.e. the stacking 
is again close-packed cubic. They still contain the ‘chemical’ stacking 
fault. This is however difficult to detect since the scattering factor for 
electrons of Sb and Te is very nearly equal, in other words electrons do 
not distinguish between Sb and Te. 
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On the same basis a plausible explanation for the particular geometry 
of the pairs of loops can be provided. As a consequence of the elastic 
interaction between vacancies and dislocations the centre of a prismatic 
loop is the preferred site for nucleation of a new prismatic loop (Frank 1956). 

It is reasonable to assume that the nucleation will take place in the 
Te, layer which is nearest to the first loop plane, i.e. in the plane indicated 
by the arrow in fig. 4(b). This loop will evidently have the same character 
as the first. The question now arises: which will be the Burgers vector 
of the partial nucleated to remove the stacking fault? As in the first 
loop, the partial has to transform an a-layer into a c-layer and it has 
therefore the choice between three possible vectors, one of which is the 
same as in the first loop. A dislocation with this vector will however not 
be favoured as compared to the others because it would be repelled by 
the segments of the existing loop. Choice of the other vectors will on 
the other hand lead to an attracting partial dislocation. The direction of 
elongation of the second loop will thus differ from the one of the first by 
the observed angle (fig. 5). 


Fig. 5 


bt 
The projection of the Burgers vectors for the two concentric loops. 


It is clear that the deviations from stoichiometry might also be accom- 
modated by the introduction of antimony interstitial loops. It is however 
well known that the geometrical stacking fault resulting from the intro- 
duction of a layer of interstitials in a cubic stacking, cannot be removed by 
glide along the plane of the loop. We should therefore find a ‘ geometrical ’ 
as well as a ‘chemical’ fault; and this should give rise to stacking-fault 
contrast, which is not observed. We can therefore exclude this possibility. 

A number of interactions between the loops has been observed. Such 
interactions have been analysed in detail in the case of zinc (Fourdeux 
et al. 1960, Fourdeaux et al., to be published) where the loops were however 
due to the condensation of point defects created during cold work. ; 

In fig. 1(b) (in the points marked A) interactions between the sessile 
loops and glissile dislocations, leading to pinning of the latter are clearly 
visible. The basal component of the Burgers vector of the loops may react 
with the complete glissile dislocation and produce a different partial 
dislocation. Using the notation introduced previously (Delavignette and 
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Amelinckx 1960) and ignoring the vertical component (¢/15) of the Burgers 
vector of the loop, we can represent the reaction as shown in fig. 6 (a). 
Along the segments X ZY the reaction 
yC St CB+yB 

takes place. 

Repulsion between the glissile dislocations and the loop may also occur 
(fig. 6(b)) in fig. 2, e.g. the dislocation UVW is clearly repelled by the 
large loop, whilst other dislocations are attracted. 


Fig. 6 
C 
C C 
B X Y B 
7é 
(a) (d) 


Interaction between a sessile loop and a glissile basal dislocation. (a) Attraction; 
(6) repulsion. 


4.3. Equilibrium Shape of the Dislocation Loops 

The isolated large loops that result after slow cooling may be assumed 
to have adopted an equilibrium shape. ° 

We will now calculate this equilibrium shape for a loop of the kind 
considered in the previous paragraphs. 

Its shape will be described by the two parameters x and h (fig. 7). The 
total area of the loop is then 

aa (Ee. ech Bee a 
We assume that the equilibrium shape is such that for constant area S, 
the total energy # is a minimum. 

We call the magnitude of the Burgers vector 6; AB and ED are pure 
edges with 6 as a vector. The edge components of the mixed segments 
are given by 

be? =(t)6*(sm*a4+4c0e"a) 5 Ps) 
where a is the angle between the normal to the loop plane and the Burgers 
vector. The screw component of the mixed segment is 


bs" =(*)b*sin?a. . . &(3) 
The line energy per unit length of the sige ane? AB dae ED i is 
He = [ub*/4er(1—y) In yr), 9. 


whilst the corresponding quantity for the mixed sections is: 
Em = (u/47) In (71/79) [be?/(1—v) +052]. 2. 2 2 7... (5) 
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Neglecting the interaction energy between different segments of the loop 
the total energy associated with the loop is given by 
H=2x He+4yHy. . . . . . . e (6) 


If we assume the stacking to be perfectly close-packed cubic we have 
sin? a= % and cos?«=+; if we further accept v=4 one finds 


H=2 jo In (3) [4mt—») | [w+ (3)y]. Meee ee ia (d.) 

ro 
Since x=(S/h)—(/3/6)h and y=(4)hy/3, one can easily minimize E 
with respect to h, and subject to the condition S=constant, this leads to 
h? =188/74/3. PCS te eee 


The quantity w/h, which is more easily compared with the observations, 
is then given by 
Mi heats 22) 20-9 Oot cat yon 8 bos aac eet (9) 


To illustrate the notation used in calculating the equilibrium shape of the loops. 


A number of well-formed loops were measured. From the deviations 
of the angles from 60°, the inclination of the plane of the foil with respect 
to the lectron beam can be calculated and hence taken into account in 
the measurments of h and w. 

The results are given in the table. The agreement is not too bad especially 
since we neglected the interaction energy ; it should be best for large loops. 
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These observations lend further support to the model adopted for the 
loops; no conclusion concerning the generation mechanism can of course 
be drawn from these considerations. 


wt h(t) wht 

22 27 a 

27 35 0-77 

21 25 0-84 fig. 1 (¢) 
22 o7 0-81 

18 23 0-78 

42 52 0-81 


+ In arbitrary units, made on an enlarged photograph. 
t In the regular hexagon w/h=1-16. 


§ 5. CONCLUSIONS 


It is clear that the explanation offered here is somewhat speculative 
in its details. The important point which we hope to have made plausible 
is that the deviations from stoichiometry can be accommodated in a 
crystal, at least partly, by the precipitation of a number of vacancies in 
certain lattice planes into dislocation loops. Especially in layer structures 
this should be a probable feature. This picture may account for the large 
influence of small changes in stoichiometric composition on the mechanical 
properties of some of these compounds. The sessile loops should indeed 
very much hamper the movement of dislocation on the basal plane. 
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The Scattering of Phonons by Bound Electrons in a Semiconductor} 


By Gy Pyin 
Cavendish Laboratory, University of Cambridget 


[Received November 28, 1960] 


ABSTRACT 


In a fairly pure semiconductor at low temperatures, the carriers are in states 
bound to impurity atoms. In compensated material, there will be some 
bound states which are not occupied. 

A carrier can jump from one impurity state to an unoccupied neighbour, 
and thereby interact with lattice waves, setting up a thermal resistance. We 
calculate the effective mean free path of phonons affected by this mechanism. 
This mean free path varies inversely as the number of carriers and the number 
of vacancies, and as a function of temperature has a minimum, at about 0-:2°x 
for ‘light’ holes in germanium. Some phenomena in the thermal con- 
ductivity of germanium are thought to be due to this effect. 


THE thermal conductivity of germanium has been measured at low 
_temperatures by Carruthers et al. (1957), and Carruthers (private com- 
munication, 1958). The measurements were made at temperatures 
between 0-2 and 90°K, on a number of specimens of different purities, which 
contained from 101° to 101° p-type carriers at room temperature. 

Under these conditions, the heat is carried almost entirely by lattice 
waves (Rosenberg 1954), and in order to explain the results we must 
‘consider the various mechanisms by which phonons may be scattered. 
Point defects in the lattice (isotopes and chemical impurities) cause 
Rayleigh scattering, which can be estimated by extrapolation from the 
high-temperature results. The effect of boundary scattering can also be 
calculated (Casimir 1938), and approximately allowed for. The residual 
thermal resistance must be due to the scattering of phonons by some other 
mechanism. This resistance was found to depend strongly on the impurity 
content of the specimen, and is given in fig. 1 expressed as an effective 
mean free path for phonons. 

In the purest specimens (Ge 2, 3, 4, 5 and 10, containing up to 1015 
impurity atoms per cm’), phonons have relatively long mean free paths, 
which indicate that the thermal resistance is fairly well explained by 
‘boundary and isotope scattering. But for the less pure specimens there 
is a significant amount of scattering of phonons, which indicates the 


presence of a new mechanism. 


+ Communicated by Dr. J.M. Ziman. 
t Present address: Theoretical Physics Division, A-E.R.E., Harwell. 
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The behaviour of the most impure specimens (Ge 11 and 12, which 
contained more than 10! impurity atoms per cm?) has been explained 
by Ziman (1956b). The additional thermal Tesistanee is attributed to the 
absorption and emission of phonons by ‘free’ electrons which behave 
like a degenerate Fermi gas because of the strong overlap ees centres. 

One specimen, of intermediate purity (Ge 7, with 2-3 x 101 impurity ce 
per cm3), behaved in a very surprising manner, quite different from either 0 
the extreme types. The distinctive feature of the results for this specimen is 
that the effective mean free path decreases as the temperature falls, from RO: 
5°K down to the lowest temperatures at which observations were made (0-2°R). 
This contrasts markedly with the ‘normal’ increase, which would be caused 
by Rayleigh scattering, and also with the type of curve predicted for va 
impure specimens, which has a minimum mean free path at about 5°K. 


Fig. 1 


DENSITY OF MEAN FREE 
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*® FOR Ge2 THE EFFECT OF 
BOUNDARY SCATTERING IS NOT 


Ge2 ELIMINATED. 


0-2 ie) ite) 
TCK) 

The positive value of 0A/07' implies that the interaction increases as the 
phonon wavelength increases, ie. as the local distortion of the lattice 
becomes smaller. It seems unlikely that there will be a continuing increase 
for even longer wavelength phonons, and we therefore presume that a 
minimum exists in the curve, at or below 0-2°K; the scatter of the experi- 


menta] points is very large at the lowest temperatures, and it is extremely 
difficult to detect a minimum. 


Si de t 


Scattering of Phonons by Bound Electrons in a Semiconductor 611 


In the present paper a calculation is made to try to explain this anomaly. 
The additional thermal resistance is still assumed to be due to interaction 
with electrons introduced into the semiconductor by impurities, but 
because of the smaller concentration of impurities, one no longer has a 
degenerate gas, but electrons are able to tunnel from one centre to another. 
This is only possible if the semiconductor is compensated (Conwell 1956) 
so that there are vacant centres for the electron to tunnel into. As pointed 
out by Mott (1956) the charged minority carries introduce a considerable 
random field, so that the potential energy of some centres is lower than 
others. An activation energy is thus necessary to allow the electrons to 
move, and this accounts for the temperature-dependence of the resistivity. 
We shall, nonetheless, assume that there is a large number of pairs of 
centres which have nearly the same energy, andwhich share a single electron, 
so that resonance between them is possible. This may be possible, even 
at quite low temperatures (<1°K) if the degree of compensation is high. 
We do not know that what degree of compensation was in the experiments 
quoted ; further experimental work in this would be of interest. We shall 
assume that acceptors and donors are present in comparable concentration. 

Our model is thus that of a number of pairs of centres distant R apart, 
each pair sharing an electron. The energy separation of the two lowest 
states, similar to the bonding and anti-bonding states of H,+, depends 
sensitively on R. The phonon is scattered due to a process in which the 
electron jumps from one of these states to the other. We calculate the 
probability of this process; the actual distribution of ions in the crystal 
is taken into account by averaging this result over all values of R. 

In order to obtain this average, we need to know the distribution of 
separations between pairs of ions. The most probable separation A, in 
material containing V impurity atoms per cm? is given by 


AN Rei Laaeu ae airs Md et as dae 1) 


We define the relative separation r=R/R,, the probability distribution 


for which is 
ta EO Re wel ee OE BS 


The wave functions for the simple system of two ions and one electron 
are analogous to these of a hydrogen molecule ion H,* (Pauling and Wilson 
1935). If we assume that the basic lattice is a continuum with dielectric 
constant ¢ (of value 16 for germanium), and that the electrons have 
effective mass m*, then the units of length and energy corresponding to 
the Bohr radius a, and energy W, are 


DAPSCUDN eee Oe a C) 
Wie pane ey clipes. sere, (4) 


Since the temperature is low, we consider only the two lowest eigenstates 
of the system. These are the symmetrical and antisymmetrical combina- 


tions w,(r) of the hydrogen-like wave functions u(r + R/2) for the ground 
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states of the individual ions. These eigenstates have energies £,, given by 


1 ce el, 
a Wea 23 Me pa eaten arenes iste pe 
PEA b+ a 1+A \ 
where J=1-—(14+ R)exp(—28£), 


K=(14+ R)exp(—8), 

A=(14+ 2+ R?/3) exp (— R), 
the units for R and E being a and W respectively. The difference between 
the energies of the two states is 


E(R) =: 


) 
he cae ie 


We write 

E..=E mean F $1 (£). 

The wave function for the unperturbed system consists of a linear 
combination of the eigenfunctions w., with their respective energies H , 
[g..w, exp {iH(R)t/2h} +g _w_ exp {—tH(R)t/2h}] 

x exp {—iE meant /h}. 
This represents the quantum mechanical resonance of the electron between 
the two impurity atoms, with angular frequency H(R)/2h. The probability 
that an electron is in one of the states w, depends on its energy, and is 
given by the Fermi—Dirac distribution : 


lgx[PP=fe=([1+ exp (F/eT)}>: 


The effect of lattice waves on the system is to cause transitions between 
the two eigenstates, with the simultaneous absorption and emission 
of phonons. The present situtation is rather unusual, because two eigen- 
states have non-zero initial probability amplitudes: 


a(t=0)g,exp(ip), a_(t=0)g_exp (—14) 
where ¢ is an arbitrary phase. Conventional perturabtion theory results 
are based on the assumption that the system is initially completely in a 
definite eigenstate. The modification to cover the present case needs a 
little care, but is not difficult. It will be noted that since the squares of 
the probability amplitudes are involved, the arbitrary phase ¢ in the initial 
probability amplitude cancels out. 

Using first-order perturbation theory, and the deformation potential 
method, with coupling constant G=2-3ev (Shockley 1950), we obtain 
expressions for the probability amplitudes a _,(¢) of the eigenstates at time ¢, 
from which are derived the transition probabilities. The result of this 
calculation is that the process of emission of a phonon p~, together with 


an electronic transition from the antisymmetric to the symmetric state 
occurs at the me 


Weep: —) =f.) SED (4 p— yP3H(R) ito), (1) 


where mp and wy are the number and frequency of phonons in mode Pr. 
The crystal, of volume V, is assumed to contain N atoms of mass M: 
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thus its density is p=MN/V. The matrix element (+p— » depends on. 
the wave functions of the electronic states : 


{(+p—)= [ero exp (—ip.r)w_(r)dr 


=i[1+(4ap)?]-* sin p. R. Seen ee ee sate 5 (8) 
This result for the interaction at a single pair of impurity atoms is now 
averaged over all values of the inter-ion separation R. Since each value 
of R corresponds to a particular value of E(R), the delta function is 
removed, and instead we have a relation between R and p which gives 
the energy conservation condition for the interaction : 
E(R) =fhwp. eNO ote per ies cecihl (0) 
Thus pairs of ions which are distant from each other are responsible for 
scattering low-energy phonons, and close pairs of ions scatter high-energy 
phonons. 
The total effective mean free path / is determined by the method of 
Ziman (1956 a), which is based on the assumption that N-processes are 
dominant in determining the phonon distribution. We write 


1=19,-A,(T), 
where 
1 W pv? 
eer . Sea . . . . . . 10 
fon n(N—n) ° 4n*G?a?’ an 
and 
Pee are? ee — 1h Sap 
A, 1674 if pe ee) tap) 2 eR) 
ee a eR RRs) dan oe cogs Oude TT) 


in which «=hw,/kT, E(R)=hw, =hvp, E'(R)=dk/dk. 

A graph of the function A,(7') (obtained by numerical integration) is 
given in fig. 2 (curve 1). The essential form of the curve is due to the 
differing interactions of the dominant phonons at different temperatures. 

At high temperatures, the phonon wavelength is short compared with 
the size of the wave function for a bound electron on a single impurity, 
and adjacent half-wavelengths of the phonon give contributions which 
cancel in the matrix element (8). On the other hand, at low temperatures 
the dominant phonon wavelength is long; thus the basic lattice is only 
slightly distorted, causing a relatively small de.ormation potential, and 
little interaction with the electron. The interaction is greatest when the 
dominant phonon wavelength is somewhat greater than the extent of the 
electronic wave function for a single ion. 

The density of impurity ions does not strongly affect this result, because 
the distribution (2) of inter-ion distances does not change with impurity 
density. Of course, the actual magnitude of the interaction does depend 
on the impurity density. 

We now consider how well the above calculation agrees with the observed 
results for specimen Ge 7. This involves putting numerical values for the 


614 I. C. Pyle on the 


various quantities into the expression (10) and finding whether the results 
are consistent with a mean free path of the observed magnitude decreasing 
as the temperature falls to 0:2°K. 


Fig. 2 


irexe) 


I 10 100 


The temperature at which the minimum mean free path is predicted 
depends on the effective mass parameter m*/m, and not on the density 
of impurities VN. But the actual situation in germanium is not as simple 
as has been assumed for this calculation. The effective mass of holes in 
germanium is not a definite, well defined quantity, but has a more com- 
plicated structure than the simple scalar here assumed. Due to degeneracy, 
there are two types of holes in germanium: ‘ heavy’ holes, with effective 
mass 0-3m, and ‘light’ holes, for which m*=0-04m (Dresselhaus et al. 
1953). If we assume that heavy and light holes interact independently 
with the phonons, there should be minima in the mean free path curve 
at temperatures of 1-5 and 0-2°K. 

Both these values are lower than the temperature of the minimum 
calculated by Ziman (1956 b), for the impurity band model. The value for 
light holes is consistent with the observed results. 

The magnitude of the minimum mean free path depends on the quantity 
lo, which involves both the effective mass and the impurity density. 
Substituting numerical values for the other quantities, we obtain 


1 1 m* 3 36 
n= a(S) 1:36 <.10°* cm, eo een ee 


The density of impurity atoms N, is not known accurately, because the 
only measurements on the specimens indicating their purity was their room 
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temperature electrical conductivity, which gives the number of holes which 
are not trapped by acceptorions. However, these quantities are of thesame 
order of magnitude, and we obtain an estimate of the required result by 
assuming that there are as many vacant impurity states as there are occupied 
ones. This assumption leads to the following calculated values for the minimum 
mean free path in a sample having the composition of specimen Ge 7: 


for heavy holes 10 cm, 
for light holes 0-08 cm. 


Thus this mechanism alone would produce a phonon mean free path 
with two minima, a low one at a mean free path of 0-08 cm, and a tempera- 
ture of 1:5°K. Since the mechanism does not act in isolation, the separa- 
tion of the mean free path curve into components for different mechanisms 
is only approximate, and not reliable for large mean free paths, it is to 
be expected that only the lower minimum should be observable, and that 
the other should be masked by the general background of interactions. 

We now have an explanation of the continued fall of the mean free path 
curve down to 0-2°K, but the calculated mean free path is about ten times 
larger than was the observed. Further, the gradients of the curves are not 
the same: the observed gradient is appreciably smaller in magnitude than 
the theoretical value for temperatures about one decade above the minimum. 

We therefore consider whether further terms in the perturbation expan- 
sion lead to a larger result. Using second-order perturbation theory, we 
calculate the mean free path for phonons which are scattered at an impurity 
centre of the type described. In this case, a phonon is virtually absorbed 
by an electron, and then re-emitted in another direction. Since energy 
need not be conserved in virtual transitions, there is no condition corres- 
ponding to (9) for this interaction: all phonons may be scattered at any 
impurity centre, instead of only those which have the right energy. 

The calculation is very similar to that described above, with the same matrix 
elements, now taken in pairs, for the two virtual transitions involved in the 
scattering. The evaluation of the integrals is straightforward, without the 
rather awkward relation (9) between variables. The result can be written 


l= Ip. .A,(7’) 


where es ee Mee te eae (la) 
n(N—n) 87Ga 
and A,(7') is a complicated function, the graph of which is given in fig. 2 
(curve 2). This function is very similar to that obtained by first-order 
perturbation theory, and consequently only the magnitude of the effect 
| may be modified, without any serious change in its essential nature. 
To determine whether this second-order effect is appreciable, we must 
calculate the numerical value of J),. We find — 


1 
ee a i) ee OS LO om hy On (14) 
- which leads to the conclusion that for effective masses of order unity 
(m*/m> 0-2), the second-order interaction is stronger, but for the low 
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effective masses with which we are here concerned, the first-order inter- 
action is dominant (by a factor of 54 for m*/m=0-04): the real absorption 
and emission of phonons constitutes a much stronger mechanism than 
virtual transitions leading to scattering. 

The dependence of these results on the impurity concentration in the 
specimen is seen from the eqns. (12) and (14) which define J). For small 
concentrations of impurities, there is little interaction and a large mean 
free path. The mean free path falls as the concentration increases, as 
- long as the model is valid. But of course, for concentrations over a certain 
value, the electronic wave functions are no longer localized, and the 
electrons form a degenerate gas. This model suggests that the greatest 
thermal resistance of the described type will be shown by the specimens 
whose impurity content is just insufficient for this to occur. 

The conclusion we reach is that the model described goes some way 
towards explaining the apparent anomaly in the residual thermal resistances 
in germanium, but fails to produce detailed agreement. 

Although the present work was carried out with a particular physical 
system in mind (impurity atoms in a crystal of germanium), the same 
model would apply to any crystal containing charged impurities. Whether 
electrons or holes are bound to impurity centres, the carriers diffuse through 
the crystal by the method described, and similarly interact with lattice 
waves in the crystal. 

Pohl (1960) has shown that the thermal resistance of crystals of lithium 
fluoride containing F-centres (halogen vacancies to which electrons are 
bound) is not completely explained by theories which do not take into 
account the motion of electrons between centres. Since frozen-in vacancies 
will here always lead to compensation, the present model may help in 
explaining the discrepancy. 
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ABSTRACT 


The problem of the states of electrons moving under the influence of 
disordered scattering centres is discussed in the simplest case of one dimension 
and weak interaction. It is shown how to define a density of levels n(Z) 
and an energy and momentum density p(k, #) and the variation of those 
functions from the case of completely disordered scattering centres to the 
completely ordered case is discussed. In the weak interaction case it is 
long-range order which dominates the problem and the effect of even slight 
disorder results in the appearance of states inside the energy gap of the 
ordered system. 


§ 1. INTRODUCTION 


In order to initiate a discussion of the electronic structure of disordered 
systems it is natural to consider firstly the simplest model which retains 
the essentials of the disordered state. So this paper will consider the 
electrons non-interacting with themselves except inasmuch as the Pauli 
principle applies, and interacting weakly with atoms which are represented 
by fixed predefined potential wells. The final simplification is to consider 
only a one-dimensional system, a restriction which is removed in a second 
paper. 

There are two limiting cases which exist when discussing order, those of 
complete order and complete disorder. In the case of complete order the 
periodicity of the atom lattice implies the periodicity of the wave functions 
and this theorem leads, irrespective of detailed calculation, to the existence 
of gaps in the energy spectrum. The other extreme is the state of complete 
disorder in which it is purposeless to discuss any particular wave function 
corresponding to any particular configuration and instead one must 
discuss the probability of an electron having a certain momentum and a 
certain energy. It will be shown that a mathematical formalism can be 
set up which, within the limitation of the model, will describe both of the 
limits and which can then be applied with confidence to the state of dis- 
order. Before presenting this formalism it is worth discussing quali- 
tatively the kinds of situation that could occur. Consider then a perfect 
lattice and a wave function with the periodicity of the lattice. Any such 
wave can be split into two components, one of which has its nodes at the 
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atoms and the other its antinodes. The energy associated with these two 
waves will be separated by the gap energy. It is not possible to have a 
wave of this frequency which has sometimes nodes and sometimes anti- 
nodes, it is a case of all or nothing. Now consider the completely dis- 
ordered case. A wave function now will not be periodic, so for comparison 
consider one which has the same number of nodes. Then sometimes such 
a wave function will, relative to the atoms, locally resemble the wave 


Fig. | Fig. 2 
n(E) 


R(E) 


Fig. 3 


n(E) 


fies 


functions whose nodes were at the atoms, sometimes the wave function 
whose antinodes were at the atoms, and sometimes it will be a rapidl 
varying mixture. Now as one moves away from the state of com he 
order to the state of complete disorder one must eventuall ee th 
chaotic wave function described above. - 
Working back to order, two kinds of situation can be envisaged. One 
can think of a wave function which is slightly disordered in a we that 
follows the lattice. But one can also consider that this ‘following the — 
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lattice’ will not last throughout the lattice but owing to the disorder it is 
possible to have functions which change over from being the kind with 
nodes at atoms to the kind with antinodes at atoms and thus having an 
intermediate energy. The ability of the slightly disordered system to 
sustain states of this latter type has to be investigated in the detailed cal- 
culations but it is clear that there are two ways in which an energy gap in 
the density of states can be filled upon the onset of disorder, for of course 
ultimately it must be filled when the system is highly disordered. The 
two ways are exemplified by figs. 2 and 3 which can be considered as 
derived from fig. 1, which shows the distribution of states in a one- 
dimensional perfect lattice, by the imposition of disorder. In fig. 2 the 
walls of the gap soften whilst still leaving a gap in between them, and it 
will be shown below that such a figure results from a state which has short 
range disorder whilst retaining long-range order. This is the state of 
affairs to be expected when the disorder merely distorts the wave functions 
of the electron. In fig. 3 the floor of the gap rises immediately giving a 
non-zero density of states right across the gap. This state of affairs will 
be shown to arise with the failure of long-range order and hence the 
possibility of getting more radically different wave functions. 

Thus the calculations below will suggest that in the free electron model 
a liquid metal, inasmuch as its properties may be suggested by the one- 
dimensional model, will have a density of states like fig. 3, since the failure 
of long-range order is the characteristic of the liquid state. 


§ 2. THe Basic ForRMULATION 
The Schrédinger equation for an electron moving in the presence of 
fixed centres R, which give rise to electron-centre potentials U(r—,) is 
i Ome 
5 -— Ve4 SU(r—R,)) ean 1 od) 5.1) 
In order to be able to discuss average properties however the solution 
G of the inhomogeneous equation, will be required, 


iO ieste, es Sate tery 
C5 = 5 V+ DU R,)) Grr a 

=d(r—r’)d(t—U) eyes ite Ud ug te aro) 
A solution G, to this equation is given in terms of the eigenfunction of the 
Schrodinger equation 


(r,75 t= 2S daledbtte exp | — ZBalt—1)] OF-1) 2.8) 


where Oa 4 s¢ 
= ubat 
In terms of energy by Fourier transforming G one has 
G(r, 7; B)=dbalr)bn* (1 (LE — HE, +t€)™. ene (2:4) 


The ie in the denominator is needed to define the sense in which the pole 


Z2R2 
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at H—E,, is to be taken, since there are two basic solutions of (2.2) corre- 
sponding to outgoing waves, i.e. to (2.3), and also to ingoing waves which 
use —ie in the denominator. Now since the functions ¥,, are normalized 


| ee rR) dra (HB, £se\z ee a ee 

so that since : : 
(HE, —ic)1—(E-E, tie) 1=2i8(H—-E,), . . (2.6) 
| (G(r,r; Z)-Gir,r; #)|dr=— Qrid(H—-EH,) . . (2.7) 
Dan). ieee) 


where n(Z) is the density of levels. The average density of levels is obtained 
by averaging n(Z) over the configurations of centres in the probability 
with which the various configurations occur. Denoting this averaging 
process by brackets one therefore has 


— 2m (n(H))= | [K0.) — «G_)r, ee 


Now if, as will be the case, the averaging is such as to destroy all inhomo- 
geneity in the problem the averaged Green functions G, and G_ will be 
functions of r—r’ alone 


(G(r, r°; B))— (Gr, v's E)) 


= —2nip(r—1'; E). wee 
In terms of the Fourier transform of p: 
(n(B)y= |” p(k B)dk. . eocce 4 ee 


To interpret this result physically return to (2.4) and write every y,,(7r) in 
terms of its Fourier decomposition : 


dnlr) = | o(k) exp (ikr) dk. +) ds ee eee 
Then 
(G(r, 9°53 BY) = — > (a™(k)a*™(9)) exp (tkr — tr’) 
x (H— EF, +t). Ric t ee see 
But this average must be a function of r—r’, thus the average of the a’s 
must vanish unless k=j. This means that one can write 


p(k, HB) = 2nd (ol(k)a*™(k))S(H—E,). . . . (2.14) 


Physically therefore p(k, Z) corresponds to the ‘mean square’ probability 
of finding an electron having energy H and momentum k. Even though 
the main objective is to calculate n(Z), to do so via p(k, £) will greatly 
clarify the situation. Mathematical manipulations are however much 
easier to perform on the (@) rather than on p since the (G)’s satisfy 
inhomogeneous equations and so allow simple expansions as will appear 
in detail below. The analysis now to be given follows closely that given by 
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Edwards (1958) (hereafter referred to as I). Similar treatments are given 
by Abrams and Weiss (1958), Abrikosov ed al. (1959) and Langer (1960) 


Another approach is due to Chester and Thelling (1959). The basic idea 
is that G can be averaged by considering its expansion in terms of the 


deviation of the potential ( V = }U(x—R,) ) from its mean and this series 


can be averaged without difficulty. After averaging the series may be 
justifiably summed up again. 
Write symbolically 


G,=G,— | GV —<V))Q, 


+ [arr —rmedv— CG... Se ee (215) 
where G) satisfies, in one dimension explicitly, 
fe Go he O7 ; 
— ,= ye te me 
€ ot 2m dx? eae i) Oa on 
= 6(~%—2’)d(t—t’); doe Wes week Lael) 


Gy = (27th/m(t —t’))#?O(t -t’) 
+exp E (v—2’)2(t-t’)-1 — ; ryt) | Bt ee (Oak) 


or in Fourier transform upon ¢ and x (the same symbol G, will be used), 


using units which absorb the constants 
Gy =(#— +1e+<V >). Renee rii( a 8) 


In (2.15) | G,VG, means, for example, 


| Gate: yx; t,t )Viy)Goly, x; U,U) dy dt”. 
The series can now be averaged. Use brackets to denote the average 
value 


(6..)=Gy— | OV —(V YG 
+ | GV — (VK VV)... eth 75) 
Retinihe the savies’ one has 
(B—i+ier (V)— | (VGV)+ | (VIG) + ..-)(@.)=1. (2.20) 


Now it has been shown, for example in I, that when V is weak all the odd 
averages (except (V)) are small compared with the even averages which 
precede them. ‘Thus the validity of the summation above depends on ~ 
the smallness of 

(VG VG. VGV Y—(VGV GG KVGV>)  . . . (2.21) 
and similar higher combinations. ‘That this is so was explicitly shown im 
I for the case of a random distribution of weak interactions, and as will 
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appear below it gives the standard results in the completely ordered case. 
A full discussion of this point in the one-dimensional case is given in a 
recent paper by Hisenschitz and Sah (1960). However, it is not so easily 
justified in three dimensions as in one. Since it is therefore discussed in 
detail in a following paper, which deals with the three-dimensional case, it 
will not be discussed further here. Henceforward, then, (2.20) will be 
adopted as the basis of this paper. 


§ 3. THE Description OF DISORDER 


The average (U(x—R,)U(y—,)) must be a function of x—y from the 
spatial homogeneity of the problem. This can be expressed in terms of 
the correlation function c(k) defined by 


o(k)= > (exp tk(R, = Ry), 2) os eee 
a, B 
so that if V(k) is the Fourier transform of V(x) 


> (U@—R,)Uly~R,))= | exp [ik(w—y)]V (k) V*(k)e(k) dk. 


If # is redefined by displacing it to H+<V), to avoid the repeated 
appearance of (V), then (2.20) written out in full is 


| B= tie) | (B—(j—k)? +ie)AV(j)V*(je(j) a 


+] (HE —(j—k)*+ ‘A NMOT*O) | (G.Ak;f)}=1,. (3.2) 


Since V is assumed weak and there is no question of discussing bound 
states, V can be considered a 6 function in space, so that V(k) is a constant v. 
As a convenient measure of the strength of the potential we shall take 


A=0(N Lye Aes Soe aera 


where N is the number of scattering centres and LZ the total volume. 
Then (3.1) becomes 


[#- Wa ienaine | o( (EH —(k—j)2 +ie)4dj 


Hada (6, y=. ia: Pedgeuaaee 


Now consider the structure of c(k). For a completely disordered system 
for all k40, only the terms with «= B survive the averaging 


, 1 
Slexpik(R,-Ry))= 75 Y | | exp (ik(R,—R,))dB,dR, (3.5) 


+z | aR, 
=N (k#0).- ay Scie ean 


When k=0 one gets N* instead, and this is clearly always the value 
irrespective of ordering. In passing to a large number of particles in a 
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large volume, c(k) behaves like 8(k) at k=0, and this behaviour exactly 
cancels the terms in (3.1) coming from (V)G,<V). In future this term 
and the 8(k) function in c(k) w ill not be considered further since they 
always cancel, and the random case will be taken to be simply c(k) =. 
For a completely ordered system R,—R,=(L/N)(«— 8), where L is the 
length of the system, taken to have cyclic peondary conditions. Taking 
&,, as the origin, since k has to take the values m/L where m is an integer 


N-1 
> exp (tkLB/N)={1—exp (ikL)}/(1—exp(ikL/N)} . . (3.7) 
0 
=<0 if kAnk fi owrvots wy eel Fbto.8) 
=N —it- k=nk Ts ete, Geena (3.9) 


where 27xN=LK, K being the reciprocal lattice vector. In the limit of 
L and N being large so that one can treat k as a continuous variable, eqn. 
(3.7) becomes 


Eee (ok \iee a ha) 8 BA0) 


n=-— 0 


These two limits are of course both standard results. Now consider a 
disordered system. Firstly consider each atom distributed by a function 
f(R£,— RO) about its lattice point R,© which is fixed in space. Then 


> {exp tk(R, — R,)) 


=L ¥ | f(R,—R,) f(R,—R,) exp [ik(R,—R,)] dR, dR, 


a= 


+1 S| fh, ROAR, 
=L? > | f(R,)F(B,') exp [1k(B,’ ie [ik(R,O — RK) | dk, dk, 
ae +L 2 [r@ 
ite fei > f(b e(k—nk) BAL) 


n=— oO 
The effect of this kind of disordering is to alter the strength of the 6 functions 
in c(k) and to add to it a continuous function, but it does not alter the 6 
function nature. Consider next a disordering in which each atom has a 
distribution function relative to its nearest neighbour, for the sake of 
definiteness on the left say. Then 


eee hee Gene fe wt deh (3.12) 
and «, is distributed via some f(«,). This gives 
( » exp ik(R, — R,)) 
| > exp | seEN-(a—B B)—1tk x c, | Iftende (3.13) 
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Taking R, as origin as before, and using # for ‘real part of’: 
c(k) = AS exp [(ikLN-1 + log f(k))B] 
= A{1—exp [(ikLN + log f(k))N ]} 


x {1—exp [(ikDN— + log f(k))]}}™ 
= BA —f*(k) {1 —f(k) exp ((LkN)}™. = Soa 


Fig. 4 


a 


uc 


K 


Apart from the 6 function behaviour at k=0, which as has been noted 
above is unaltered and irrelevant, the 6 functions of the ordered state are 
in this case smudged out in a way depending upon f(k), for example if 
f(€) = (27)? exp (— $e*u-*), 
then 
J (k) = (277)? exp (— p42 n-*). 
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A graphical representation of c(h) for these cases is given below using the 

height given to the 8 functions as a measure of their strength and omitting 

the 6 function at the origin. : 
Fig. 6 


e 


Hig, 7 


Since 
| (A) ai=N50) 
=NL 
the area under each of these curves is the same and one can consider fig. 5 
for the perfectly ordered case becoming partly disordered into fig. 7 then 
being finally flattened into the complete disorder fig. 4. Comparison of 
fig. 6 with fig. 7 is of interest since the effect of the disorder in 
fig. 6 upon the 6 function is to alter their strength but not their widths, 
which are altered in fig. 7. The width of the ‘6’ function is thus a measure 
of the long-range disorder present in fig. 7 but not in fig. 6, whilst the back- 
ground is a measure of the short-range disorder present in both. Since 
the case of a liquid is more reasonably described by fig. 7 than fig. 6, it 
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will be case fig. 7 which is discussed in detail below, i.e. the consequences 
of the failure of long-range order. It must be emphasized however that 
in the tight-binding approximation the effects of figs. 6 and 7 are not very 
different, so that at present the analysis below only applies in the weak 
interaction case. 


§ 4. Tue Evatuation OF p(k, £) 
The kernel in eqn. (3.1) can be split into its real and imaginary parts by 
the identity 
(HE —9? +1) 1 = P(E —j?) 1 —118(£ —9?) ee ee 
where P indicates that the singularity is understood as a Cauchy principal 
part. Then 


Pr 


M | ok —j(B—F+ ie) 
=P i o(k —)(H 92) —tad2e(k— +/E)[2,/B. . (4.2) 


=Q—-v1I say. Pon tienen lca tna ge 


In the completely disordered case c is a constant and the principal part 
(for this case of 5 function potentials) vanishes and the kernel is pure 
imaginary 


['saNA204/7Ro oS (oes ee 
Thus (G5s(-—BPHl. 2 ee 
and o(k, E)=T(#-eY4+T2}4, =.) 2 ee 


a Poisson distribution about the original parabola H = k? which tends back 
to 6(H —k?) as [ tends to zero. 

In the completely ordered system it is only the principal part which 
matters since the imaginary part is a 6 function which vanishes where p 
(k, H) is non-vanishing. The value of p (k, £) is itself a 5 function 


np(k, E) =I (z- 2 +ie—B5(E (k—nkyey2)” 


=n (BINS (B—(k—nb)). 2 aan 


Since A is small one may consider the region where only one of the denomi- 
natorg in the sum matters, »=1, giving in this neighbourhood 


i, p(k, #) =8([(E — k*)(# — (k— K)®) —2][B —(k—K)?]}). (4.8) 
Ey= (k =~ Ky 
By = 3[k* + (k—K)? + /{[k?—(K —k)?}? + 4033] (4.9) 


By =3[k? + (k—K)— s/{[k® —(K— ke + 4023] 
this can be written 
p(k, £) =0((H—H,]|#-E,)|[H-H,}>) 
= (#— E,)(E— E,)~18(E — Es) 
+(H—#,)(H-E,)8(E—-H,). . . . (4,10) 
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So the &, # diagram consists of the curve H=E, and =H, and the work 
above is reproducing the usual weak interaction theory as for example 
given by Peierls (1955). If one plots H=k? (=H, say), H=H, =, 
=H, as on fig. 8, one finds H=H, and E=E, tend to Hey, =H at 
their extremities and have the usual energy gap of A, as shown in fig. 8. 


Fig. 8 


Dy 


See , 


The probability function p(k, #) (which in this case, since it is entirely 
made up of 6 functions, is really a ‘certainty’ function) tells one more 
than just the k, # relation, for the factors (H—H#,) (f—E,)“ and 
(H—E,) (E—£,)— give the density of states along the curves =H and 
H=E£,. By adding in all the other values of 7 in the sum in (4.7) one 
obtains further energy gaps where the other parabolae H =(k—nK)? meet 
the principal parabola H=E,, and by going to approximations higher 
than (3.2), also the energy gaps where these parabolae meet one another 
until the complete set of Bloch curves is developed. The additional 
information in p(k, £) is the amount of each wave present, i.e. the magnitude 
of the fourier coefficients when a Bloch wave is analysed into plane waves 
whose frequencies of course differ by multiples of K. It is useful to 
visualize this additional information by considering the three-dimensional 
surface p=p(k, #), fig. 8 being the projection of this surface onto the #, k 
plane. The extra information contained in the values of p, e.g. in the 
perfectly ordered case the strength of the 6 function, is lost if one considers 
only a basic cell. In addition one does not wish to invoke the Bloch 
periodicity theorem since this only holds in the case of a perfectly 
periodic system. The p(k, #) diagram for the completely disordered 
system consists of a Poisson distribution based upon the original parabola 
E=E,. According to the formula (4.4), I is infinite at H=0 but the 
approximation (3.1) is not valid here since the scattering by a single 
potential well cannot be treated by Born approximation near #=0. A 
more accurate treatment using the exact scattering modifies H~'? to 
E-2(1 + 34AH-12)-1 (as will be discussed in a subsequent paper on the 
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tight binding approximation) so in drawing the projection of p(k, &) on 
the E, k plane the width will be left finite; this point has no bearing on 
the rest of the paper. The diagram is given in fig. 9. 


Fig. 9 


This diagram does not of course have a sharp edge to its width but it 
is difficult to draw otherwise. Since the decrease in order from fig. 5 to 
fig. 4 goes via fig. 7 there must be a p(k, H) diagram corresponding to fig. 7 
which lies between figs. 5 and 4. The key to understanding the situation 
can be obtained by looking at the simple function 


which is shown in fig. 10. 
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Now consider y averaged by the definition 


fie\de 
yY)j=P |= 
(yy=P | AO 
in particular consider the two examples 
Un d 
(y)= 5s | : 
2a J _, l—x+e 


1 
= 5, log (1-2 —a)/(1—x+w)], 


which is given in fig. 11 and 
a * de 1 
==. P. ie Sse ee SE 
CPi ie (2402) (l—rre) 
=(1—2)[(1—2)?+02]4 


which is given in fig. 12. 


As « tends to zero, these diagrams tend back to y=(1—~) and at «=0 
that part of the curve which links the peaks on the left and on the right 
of the singularity at 2=1 vanishes, so the diagrams can be described as 
being very similar to the y=(1—)~1 curve but with a linking curve which 
joins the two peaks. These results can immediately be translated into 


the p(k, #) diagram, for 
pase hee Oye ade ge es ca. (4.11) 
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so the projection of p onto the k, # plane will have its maximum intensity 
near E—k?+Q=0. The integral (4.2) for Q has the form of the model 
(fig. 12) above, for one can write 

(H—(K—k)*)=48 48[(./E-(K—-k)4*+(VE+(K—-®)*] (4-12) 
and, if for the moment only the neighbourhood of the first gap k= K/2 is 
Nonearred (so that the second term corresponding to the gap at k= — K/2 
can be ignored), the models (figs. 11 and 12) above correspond to 

c(k)=0, |k—K|< 2c 


= (2a); |4— KK 20 a eee 
and 
c(k)=mal(k—K)®*Poty* 2... ow (ae 
respectively 
Fig. 12 
b, 
1 =e 


The diagrams of H — k? + Q = 0 for these two models of c(k) are respectively 
figs. 14and 15. The projection of p(k, 2) wpon the (k, Z) plane is sketched, 
as was fig. 9,in fig. 15. Again of course there are not really sharp boundaries 
to the distribution of states but the intensity dies away in a Poisson distri- 
bution manner. For model (4.13) the width exists only within the two 
parabolae H=(k—K-+.«)? and H=(k—K—.«)? which are dotted in, and 
thus only the linking branch has a width as shown in fig. 16. One can 
visualize the development of fig. 8, from the parabola H =F, by thinking 
of the parabola as being made of rubber, clasped at the two points at 
which it is cut by the curves H=(K—k+o)? and B=(K-— k—«)? and 
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pulled out along these two parabolae in opposite directions. As « tends 
to zero the amount of the curve H = #, which is to be sketched out into the 
linking curve tends to zero and the number of states in the gap vanishes. 
The number of states available at say H =k?/4, the middle of the gap, can 
be estimated when « is small by integrating (2.11) and if N, is the density 


Fig. 13 


Thy 


Fig. 14 


k 


of states at this point on the unperturbed parabola H = Hy, i.e. Ny= (2 £)™, 
Eg is the energy gap, and LH, is the energy associated with the momentum 
spread a, i.e. a”, then 


Meap=NoH,/H, (H,«<Hs). 


_ Diagrams for n(£) will be discussed in the next section. As a becomes 


large the curve (fig. 15) straightens out, the linking curve becoming parallel 
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to the E axis at \~aK and eventually becoming the diffuse parabola of 
fig. 9 as a tends to infinity. This latter limit is of course unreasonable 
since the spreading of the one peak of c(k) which has been considered so 
far will meet the others to amalgamate into fig. 4, but again of course 
leading to fig. 9. 


Fig. 15 


ny 


Fig. 16 


The linking curve is a characteristic of the disorder of fig.7. An analysis 
of the disorder represented by fig. 6 will not have this feature. In this 
case the kernel of (3.2) is like a mixture of the complete disorder kernel 
weighted with (1—f(k)f(—z)) and the completely ordered case weighted 
with f(k)f(—k). So the behaviour of the (2, k) curves will resemble the 
completely ordered case (fig. 8) with A®f(k)f(—k) replacing the inter- 
action strength A’, and the curves will also acquire a width. As f tends — 
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to zero the curves H, and H, tend to H, and EZ, and the density of states 
along H, tends to zero, so that finally one reaches fig. 9. In fact the 
behaviour is just like that of a perfectly ordered system_as the interaction 
tends to zero, with the addition of a steadily increasing width. 


Fig. 17 


§5. Tae Leven Density or THE DISORDERED SYSTEM 
_ The p(k, #) diagram must be symmetric in k and this modifies figs. 8-16 
in the states of complete and near complete order in a way which will be 
- important for the final discussion of this section, so the modified 
 #E—k2+Q=0 curves are now given. Figure 17 is tor complete order. 
J 
e. 


2s 
P.M. 


634 S, F. Edwards on the 


The limiting linking curve which in this case has zero density, has been 
dotted in. Figure 18 is for almost complete order, fig. 19 is the inter- 
mediate case between fig. 18 and fig. 17 in which the upper and linking 
branches meet at a cusp. The curve has one branch in fig. 20 in which 
there is less order, whilst complete disorder returns again to fig. 9. 


Bigwld 


E 


Fig. 20 


|< 


k 


It will be seen that the main effect is that the linking curves link with 
one another rather than with the upper principal curves which themselves 
link with one another. These curves result from using both terms in 
(4.12). Note that there is also a further branch of H—k?+Q=0 lying 
above the others and touching the higher of the two where it meets the 
axis. This branch is omitted here as it has no bearing on the analysis 
and belongs to a discussion of the higher Bloch curves. The linking 
curve fig. 17 meets the # axis at the point (p) where H= K2. The minimum 
of the upper (full) curve is at (7) and the maximum of the lower at (r), 
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these points corresponding to the energies E,, E,, both of which have 
k=K/2. These gap energies satisfy 


E— K?/4~)(E — K?/4)-; ge ee sco 
: ee 
1.e. LE = kh? +4, } (5.2) 
E,,=+k?—2. 


Thus in the spirit of the almost free electron model with a small inter- 
action A, it is clear that L,,>H,>E,, as has been drawn. Now n(£) is 
defined by 


nye | p(k, H) dk. 


To visualize it from the diagrams think of the projection of p=p(k, E) 
onto the #, k plane, itself projected onto the # axis. In the completely 
ordered case, considering only the neighbourhood of k= K/2, one gets 
from (4.10) 

n(E) oc (H — 4k? —))-12 


and similarly near H = E,; so that n(#) behaves as is shown in fig. 1 of the 
introduction. The peaks are characteristic of one dimension only. By 
looking at figs. 14 and 16 one sees the effects of the disorder to be firstly 
to round off the peaks, as the width is not zero at the peaks H,, H,; but this 
rounding is mild since the width at H,, H, is small, indeed in fig. 16 the 
peaks remain infinitely sharp. Secondly, a roughly uniform density of 
states appears between #, and H,, this being the projection of the linking 
curve on the # axis, and the resulting distribution is that of fig. 3 of the 
introduction. 

The analysis of this paper suggests then, that under the conditions of 
the almost free electron model the energy gap is immediately destroyed 
upon the slightest disorder of the system. As the disorder tends to be 
complete n(#) must tend to be a smooth curve, and the current work 
suggests this comes about by the ‘floor’ of the gap rising as the peaks 
diminish. It does not bear out another way one could think of degrading 
the gap of fig. 1, by softening the walls as in fig. 2 whilst still having the 
gap until considerable disorder has set in. This latter picture can only 
be obtained from a disorder which maintains long range order, i.e. that of 
fig. 6. The results of this paper do not agree with numerical results given 
by Makinson and Roberts (1960). Presumably the density of states in 
the gap predicted here is too small to be unambiguously discernible in the 
cases they consider. 

It must be emphasized that these results all refer to the almost free 
electron model. It is well known however that the results of the tight- 
binding approximation for a perfectly periodic lattice are topologically 
identical with the weak-binding case, so it is worth while considering 
pushing the above calculations outside their region of strict validity to 
see if they change qualitatively. Consider again then the case of 
infinitesimal disorder which leads to a linking curve (fig. 17). One 
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can think of an increasing strength of interaction pushing the point 
(p) down until it lies below the point (¢). 

In this case the density of states, being the projection of these curves 
onto the E axis, will still show a gap when disorder appears, and this gap 
will close upon increasing the disorder by closing the H,, #,, gap in fig. 21 
and by the increasing width associated with the disorder. In the almost 
free electron model the width is always large since all N atoms contribute, 


Fig. 21 


E 


Fig. 22 


A 
c 
F é 

but in the spirit of the extrapolation to tight binding the width associated 
with fig. 21 should also be considered small and of finite extension. Under 
these circumstances the n(H) curve will close in from one side but not from 
the others as in fig. 22. 

The gap exists between G and F. With increasing disorder, or with a 


weakening of the interaction constant, the peak at B (which is exaggerated 
in fig. 22) moves past the peak at CO, and the picture of the earlier discussion 


n(é) 
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reappears. ‘The critical interaction comes when A=0(K?) which clearly 
corresponds to tight binding. To attempt to make these last considerations 
rigorous one must treat the problem from the point of view of the tight- 
binding approximation, and this will be done in a subsequent paper. 


§ 6. CoNCcLUSION 


This paper has discussed the problem of disorder within a very limited 
framework and serves as an introduction to the realistic situations of 
three dimensions and of taking tight binding into account. However a 
general result of principle seems to be born out which is that the curves in 
the (k, #) plane say of H—k?+Q=0, or more generally the surfaces 
p=p(k, £), are continuously distorted by the disorder and can be under- 
stood much better when considered as a whole. The point can be made 


iy 


Fig. 23 


even more strongly if one builds in the symmetries which are available in 
all degrees of order. In the one dimensional model there is only one: 
p(k, E) is even in k, so that p=p(k?, #). If one writes k*=q and studies 
p=p(q, £) then by allowing q to be negative the curves figs. 17-20 in the 
neighbourhood of #, become simpler to understand (fig. 23). 

This is also true of p itself. Since it is the imaginary part of a function 
(@) which is analytic except for a certain number of poles and branch 
cuts, it is simpler to consider the complex (CG), q, # spaces. In three 
dimensions this geometric thinking will become very difficult, but the 
advantage of considering the distortion of surfaces which remain topo- 
logically invariant except in limiting cases is very great, and it is hoped to 
bear this point out in subsequent work. 
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ABSTRACT 


In the present paper we collect the available evidence from work-hardening 
curves, slip-line observations, and ferromagnetic measurements on the easy 
glide region of f.c.c. metals and alloys and the low-temperature work-hardening 
of hexagonal metals. New experimental results are reported on copper single 
erystals deformed at liquid oxygen temperature. A dislocation theory of 
work-hardening is given which accounts rather well for the experimental 
facts. The work-hardening rate can be calculated from slip-line data 
(distances between and lengths of slip-lines) and is found to agree well with 
the measurements. The theory is a statistical one and is based on the fact that 
in the deformation stage considered here the important stress-fields are those 
of individual dislocations rather than those of piled-up groups. The present 
paper is another example for the usefulness of slip-line studies in theories of 
work-hardening. 


§ 1. InTRODUCTION 


In a previous paper (Seeger et al. 1957) a general picture was developed 
of the dislocation processes responsible for the principal features of 
work-hardening of f.c.c. metal single crystals. A quantitative theory 
was given for the rate of rapid work-hardening in the so-called stage II 
of the stress-strain curve. In the meantime this theory was confirmed 
quantitatively by electron microscope observations of active slip lines 
(Mader 1957, Kronmiiller 1959) and by ferromagnetic measurements 
(Kronmiiller 1959, Seeger and Kronmiiller 1960, Kronmiiller and Seeger 
1961). Furthermore, a quantitative interpretation was given for the onset 
of the dynamic recovery region in the stress-strain curve of f.c.c. crystals, 
i.e. the shear stress 7,;; at the beginning of stage III (Seeger ef al. 1959, 
Wolf 1960). This theory was quantitatively confirmed by experiments 
on copper, gold (Berner 1960), and nickel single crystals (Leitz 1960). 
Qualitative or semi-quantitative interpretations are available for the 
work-hardening and work-softening characteristics in stage III, for the 
transition from the easy glide region (stage [) to stage II, and for the 
origin of such structural changes as kink bands or bands of secondary 
slip (Seeger 1958, Mader and Seeger 1960, Seeger and Mader 1960). 
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What is lacking in the present picture of work-hardening of f.c.c. metals 
is a quantitative understanding of the easy glide region (Andrade and 
Henderson 1951). It was pointed out some time ago (Mott 1952) that 
the easy-glide phenomenon in f.c.c. crystals appears to be related to the 
work-hardening of hexagonal metal crystals such as zine or cadmium. 
In our previous paper (Seeger et al. 1957) we therefore postponed the 
detailed discussion of stage I until a more thorough investigation of the 
work-hardening characteristics of hexagonal metals deforming mainly 
by basal glide became available. In the meantime the plastic deformation 
of zinc single crystals was studied experimentally (Traéuble 1958, Seeger 
and Trauble 1960), using the same approach and as far as possible the 
same experimental techniques as employed in the study of f.c.c. metals. 
The combination of quantitative slip-line investigations by optical and 
electron microscopy with measurements of stress-strain curves showed 
that two principal manifestations of the work-hardening of zinc single 
crystals deforming by basal glide have to be distinguished: at low 
temperatures (<250°K) the deformation process is indeed similar to the 
easy glide of f.c.c. crystals. The slip-line pattern is essentially that of 
fine slip (Mader 1957) with slip-line lengths of the order of magnitude of 
1mm, i.e. comparable with the diameter of the specimens. At higher 
temperatures, the slip-line pattern is rather different, the slip-lines being 
grouped into rather short bands, showing the phenomenon of ‘fanning’. 
The governing processes in the high-temperature region appeared to be 
the thermally activated climb of edge dislocations and the condensation 
of point defects to form obstacles to the gliding motion of the basal 
dislocations. 

The purpose of the present paper is to close the above-mentioned gap 
in the theory of work-hardening by giving a quantitative explanation 
of the rate of work-hardening in easy glide and in the low-temperature 
region of hexagonal metals. In the present paper we shall not consider 
the high-temperature work-hardening of hexagonal metals, in which 
recovery processes are important; we shall confine ourselves to those 
cases in which the hexagonal crystals deform mainly by basal glide. 
In §2 we shall mention first some of the experimental difficulties with 
which quantitative studies of the present topic are beset, and then 
summarize the main experimental evidence available. In §3 we shall 
outline the qualitative features of the theory, stressing in particular the 
differences to the theory obtaining in stage II. Section 4 contains the 
quantitative formulation of the theory. It involves a number of 
statistical calculations, the details of which are deferred to an appendix. 


§ 2. SUMMARY OF THE EXPERIMENTAL EVIDENCE 
2.1. Basic Facts 


; The outstanding experimental observation on the easy -glide phenomenon 
is the rather small work-hardening rate. It is an order of magnitude 
smaller than in stage IT of the work-hardening curve of f.c.c. metals. 
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The work-hardening rate in the easy glide region, which we shall denote 
by @,=(dz/de);, is rather dependent on the orientation of the tensile 
axis (Diehl 1956a). We shall fix our attention on crystal orientations 
in the middle region of the stereographic standard triangle, where 6, is 
smallest and of the order 2 x 10~4G (@ shear modulus). A brief discussion 
of the orientation dependence of 6; (Seeger 1957) will be given in §5. 

For hexagonal crystals deforming mainly on one basal glide system the 
low-temperature work-hardening rate 6/G is of the same magnitude as 
the lowest values obtained for f.c.c. crystals or even somewhat smaller. 
Whereas the easy-glide work-hardening curve in f.c.c. crystals is linear 
within experimental accuracy, this is not so for the low-temperature 
work-hardening curve of zinc or cadmium crystals. However, even in 
these cases the slope of the stress-strain curve is practically constant 
over strain intervals which are considerably larger than the typical 
extension of the easy-glide region of f.c.c. crystals. Since furthermore 
the deviation from the idealized behaviour seems to depend on the way 
the crystals are prepared (see e.g. Liicke et al. 1955, and Seeger and Trauble 
1960), we shall not discuss here the strain dependence of the low temperature 
6 (see however § 5). 

A similar idealization is necessary for f.c.c. crystals. Diehl (1956 b) 
for Cu crystals and Meissner (1959) for Ni-Co crystals showed that in 
easy glide the deformation is not uniform along the specimen. The ‘local’ 
work-hardening rate may deviate up to 30% from the average taken 
over the length of the crystal. 


2.2. Particular Problems 


The experimental investigation of the low work-hardening stage has 
to face a number of difficulties which are absent in stage II. We list some 
of them here in order to see the natural limitations to the possible agreement 
between theory and experiments. 

(a) In the majority of pure f.c.c. metals deformed at room temperature 
the easy-glide region extends only over a few per cent shear strain. For 
some techniques this is not enough to do reliable experiments. We have 
remedied this by doing experiments also on a hexagonal metal (zine) and 
on f.c.c¢. alloy crystals, in particular nickel-cobalt crystals. It was shown 
that the rather extended stage of initial low work-hardening of these alloys 
is indeed closely related to the easy glide of pure metals (Meissner 1959), 
and is positively distinct from the yield elongation observed e.g. in a-brass 
single crystals (Piercy et al. 1955), which is associated with the propagation 
of Liiders bands. To some extent it is possible to get longer easy glide 
regions in pure f.c.c. metals by working at liquid air temperatures. 

(b) For both f.c.c. crystals and hexagonal crystals the slip-lines observed 
in the electron microscope are so long that it is rarely possible to measure 
- directly the length of a slip-line. It is therefore difficult to decide whether 
the active length of slip-lines (i.e. the length of slip-lies active in a small 
strain interval) varies with strain or not. Furthermore, for f.c.c. metals 
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and also for hexagonal metals at small strains, the height of the slip-lines 
is rather small, corresponding to n= 10-20 dislocations in a glide zone. 
These small values of n make accurate determinations difficult, particularly 
if one wishes to decide whether an individual slip-line grows in height 
with increasing strain or not. 

(c) Zinc, cadmium, and to a lesser extent also magnesium crystals, 
recover at room temperature. In our investigations of the low-tem- 
perature behaviour it was therefore inadmissible to interrupt the low- 
temperature tensile tests for an extended period of time in order to obtain . 
a replica or to repolish the surface, since this would have brought in 
undesirable changes in the dislocated structures of the deformed crystals. 
We could therefore use neither selected area techniques (‘ Zielpraparation ’) 


Fig. 1 
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Number » of dislocations per slip-line as a function of shear strain for zinc ™ 
crystals deformed at 90°xK. 


nor determine the active slip-line lengths and slip-line heights as a function 
of strain. We did use the selected area method, however, in order to 
study the effect of short recovery intervals on the low-temperature slip- 
line pattern, Slip-line heights as a function of strain (fig. 1) were obtained 
by investigating crystals of the same orientation that had been strained 
(without interruption) by different amounts. For details of these investi- 
gations see Seeger and Triiuble (1960). 

(d) Whereas the work-hardening rate in stage II is very reproducible 
and hardly dependent on the impurity content, the perfection of the crystal 
or the handling and previous history of the specimen, this is not so in the 
work-hardening stage with which the present paper is concerned. 
Examples for this have been given in § 2.1, 
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The preceding discussion shows that the agreement between theory 
and experiment, which we can hope to obtain, will be seriously limited by a 
number of factors. At present it appears that due to these limitations it 
will not be possible to get agreement better than within a factor of two. 

The evidence on the temperature dependence of 6,/G is not very clear 
cut. Cu single crystals (Berner 1957, Diehl and Berner 1960) show a 
decrease by almost a factor of two when going from room temperature to 
liquid nitrogen temperature, whereas Ni single crystals (Leitz 1960) hardly 
show any temperature dependence between 90°K and 370°K. A similar 
situation obtains for hexagonal crystals at low temperatures, where 6/G 
in some instances is independent of temperature and strain rate, whereas 
in others there is a slight temperature variation (see e.g. Seeger 1958). 
Since we do not aim at a better agreement between experiment and theory 
than within a factor of two, we shall leave the question of the temperature 
variation of @/G out of our considerations. 


2.3. The Experimental Approach used in the Present Investigation 


In the following sections we shall summarize our main experimental 
results. Sections (a) to (d) are based: on slip-line studies by electron 
microscopy. Sections (e) and (f) deal with flow stress and -creep-rate 
measurements and with measurements of the approach to ferromagnetic 
saturation in Ni-Co crystals. The interpretation of the latter group of 
data requires a certain amount of dislocation theory. They are, however, 
particularly useful as guides to the type of dislocation model which is 
appropriate in the easy glide region. 

(a) Due to the difficulties mentioned in § 2.2 (c) it is not possible to follow 
the history of an individual slip-line on zinc as a function of strain without 
involving recovery. We have therefore determined the slip-line heights 
after different shear strains « of zinc single crystals at 90°xK (Trauble 1958, 
Seeger and Triéuble 1960). Figure 1 shows the average number n of 
dislocations per line as a function of strain. It will be seen that n varies 
linearly with ¢ and that it extrapolates to a finite value of n) 18 at zero 
strain. . 

(b) On Zn single crystals deformed at 90°K no ending slip-lines could be 
found within the field of view of the electron microscope. This means 
that their length was comparable with the diameter of the crystal. Ona 
copper crystal deformed at room temperature, an individual slip-line 
could be followed over 0:-6mm (Mader 1957). <A statistical determination 
of the slip-line length L by counting the frequency of ending lines was 
possible on Cu and on NiCo crystals with 20% Co (Kronmiiller 1959). In 
a NiCo crystal deformed at 20°c the easy-glide region extended to «= 0-25. 

In four strain intervals in the easy-glide region we found, independently of 
the total strain, L= 1mm. 

(c) The slip plane distance ~ is defined as the average separation between 
neighbouring active slip planes (measured perpendicular to the slip planes). 
In the electron microscope x can be determined by either measuring 
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distances between slip-lines or by counting the number of slip-lines per unit 
length. On zinc single crystals deformed at 90°K it was found (‘Trauble 
1958, Seeger and Triiuble 1960) that the average slip-line distance deter- 
mined by either one of these two methods was independent of strain for 
€>0-1. Kronmiiller (1959) observed that in the easy glide region of a 
Ni-20°% Co crystal deformed at 20°c the slip plane distance x did not 
depend on the length of the strain interval and hence should be independent 
of the total strain (for typical slip-line patterns see figs. 2 and 37). 

(d) Selected area electron microscopy was carried out on copper single 
crystals deformed at liquid oxygen temperature as shown in figs. 4 and 5 
(crystallographic orientation C/14—see Diehl 1956 a). After a certain 
amount of strain the average separation x of neighbouring slip-lines was 
determined by counting the number of slip-lines per unit length. After 
additional strain increments we determined the number of additional 


Quantitative data from slip-line pattern of crystal fig. 4 


6, (kg/mm?) 


0-0688 386 600 0-88 


0-150 330 620 0-77 


0-72 


lines by comparing the electron micrographs from the same area, but 
different strains. In this way we obtained 2 as a function of strain as 
reported in the table. The lengths LZ of the slip-lines were found by 
counting the number of ending slip-lines on each of the electron micro- 
graphs. The values of 6; are calculated from x and L=L, (eqn. 4.14) 
for the purpose of checking the theory of §4. A modified and possibly 
more accurate method for determining L was used for a second crystal. 
On an electron micrograph (extending over 250cem in the direction of 
the slip-lines, magnification x 10000) we compared directly individual 
slip-lines in areas separated by a distance of 0-25mm parallel to the slip- 
line direction and counted the number of slip-lines that were visible in only 
one of the areas. After total strains of ¢;= 0-048 and «,=0-075 we found 
L=700 and «=380A, independent of strain. These observations indi- 
cate that L and # are indeed constant, provided the strain € is not too small. 

(e) Berner (1957), Diehl and Berner (1960) investigated experimentally 
how the temperature independent part 7, and the temperature dependent 
part Tg of the total flow stress 7 vary as a function of 7 or ¢ in the easy 
glide region. On copper single crystals it was found by flow-stress 
measurements involving various temperature changes in the temperature 


+ Figures 2, 3, 4 and 5 are shown as plates. 
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interval between 90°K and 295°K that 0,= dr,/de<0,, =dr,/de and also 
that 0,/@,<7s/7,. This means that the density of the dislocation forest 
threading the principal glide plane does not vary significantly during 
easy glide. This is in complete agreement with measurements of transient 
creep of aluminium and copper single crystals, which gave good agreement 
with a theory based on the assumption that the forest density is constant 
during easy glide (Seeger 1954, Michelitsch 1959). 

(f) It was shown elsewhere (Seeger and Kronmiiller 1960, Kronmiiller 
and Seeger 1961) that in ferromagnetic crystals measurements of the 
differential magnetic susceptibility y in the range of the approach to ferro- 
magnetic saturation provide quantitative information on the density 
and the arrangement of dislocations. In particular, it can be decided 
whether the internal stresses in the material are mainly due to piled-up 
dislocation groups or to isolated dislocations. The separation beyond 
which dislocations act as individuals rather than members of a group, 
depends on the magnetic field H and can therefore be varied in the experi- 
ments. 

In the easy glide region the dislocation density is rather small. In 
order to obtain quantitative measurements, it is imperative to minimize 
all the other contributions to y in the saturation range, in particular that 
of the crystalline anisotropy. In alloys of 80% Ni and 20% Co the 
crystalline anisotropy is so small that for field strengths H 2 500 Oe y is 
determined mainly by the magnetostrictive action of the dislocation stress 
fields. We have therefore chosen single crystals of that alloy for a quanti- 
tative investigation (Kronmiiller 1959). Details are given in Appendix A 
of the present paper. The result is that the correct relation between y 
and 7 is obtained, if we assume that the dislocations are distributed in the 
various glide planes as suggested by the slip-line pictures, and that in the 
glide planes the dislocations are so far apart that they act as individuals. 


§ 3. OUTLINE OF THE THEORY 


We shall at first discuss some implications of the experimental evidence 
reported in §2.3(e) and §2.3(f). The flow-stress measurements show 
that the contribution of the dislocation forest to the work-hardening rate 
in stage I is negligibly small. 6/G@ has to be accounted for by the action 
of the long-range stress fields of the dislocations. The existence of well- 
defined slip-lines shows that the dislocations are generated from their 
sources in groups of at least ten to twenty. However, if one calculates 
the internal stress as produced by such groups, one finds on quite general 
arguments similar to those outlined by Kronmiiller and Seeger (1961) 
that the work-hardening rate should be approximately that observed in 
stage II. A similar discrepancy in order of magnitude is obtained for the 
y-7 relation. We must conclude, as already indicated by the remarks at 
the end of § 2.3(f), that in the easy glide region and also in the hexagonal 
crystals the stress fields of the dislocations must be considered individually 


and not as that of piled-up groups. 
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The main task of the present section is to understand how the difference 
from stage II comes about, where the dislocations piled up towards the 
end of a slip-line act as a group. ‘The flow stress is determined by certain 
averages over the stress fields of the dislocations, and they will necessarily 
be different for stage I and for stage II. In the present section we shall 
confine ourselves to a qualitative discussion. A quantitative treatment 
will be given for stage I in the next section, whereas for stage II it may be 
found elsewhere (Seeger 1958, Kronmiiller and Seeger 1961). 

In stage II the piled-up groups of n dislocations may be treated as 
giant dislocations of strength nb. Relaxing this simplifying assumption 
(Kronmiiller and Seeger 1961) has virtually no effect on the final result. 
The reason for this is that in stage II the extension of a piled-up group is 
smaller than the distance between groups. In stage I, however, the 
opposite is true. It is therefore essential to allow for the spacial spread of 
the piled-up groups. The quantitative results of the present paper show 
in fact that it is a very good approximation to calculate the stress fields of 
the dislocations as those of isolated dislocations. 

In the subsequent discussion we use explicitly the picture that the dis- 
location multiplication during plastic deformation occurs by the Frank- 
Read mechanism. In stage II the theory is based on the following 
criterion for the flow stress: in the area of the glide plane swept out by 
the first ring emitted from a Frank—Read source before the next ring is 
formed the opposing stress (due to internal stresses or other obstacles) 
must not be greater than the applied stress. This leads to the estimate 
that the flow stress is approximately equal to the internal stress halfway 
between piled-up groups of the same sign. 

The preceding condition is obviously necessary but may not be sufficient. 
The generation of one or two dislocation rings from a number of Frank— 
Read sources would not give a macroscopic strain (implied in any flow- 
stress definition) if the rings would get held up after having travelled only 
a short distance from their source. This is in fact likely to occur, since in 
general the rings will have fo pass rather close to other dislocations. How- 
ever, by that time more dislocations may have been produced by the source 
and the corresponding stress magnification acting on the leading dis- 
location may push this dislocation past the obstacle until it gets eventually 
stuck at an unsurmountable obstacle like the head of a piled-up group or 
a Lomer—Cottrell barrier. This mechanism works well in stage II, since 
the slip distances are short anyway, and since a piled-up group can develop 
in a rather small area behind the leading dislocation of a ring system. 

In stage I, however, the preceding consideration leads to an additional 
necessary condition for the flow stress. If a gliding dislocation gets 
stuck because it has to pass too close to another dislocation, it has to wait 
until enough dislocation multiplication at the source with the accom- 
panying stress multiplication has occurred in order to push it over the 
obstacle. In fact, in stage I this appears to be the determining mechanism. 
The analysis to be given in § 4 provides an equation relating the flow stress 
to the number of dislocations generated from a source. 
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The dislocation model which we shall adopt for the quantitative treat- 
ment is the following: dislocations are generated by Frank—Read sources. 
They get finally held up by obstacles which are grown into the er ystals 
(e.g. sub-boundaries) and—as an idealization—do not increase in number 
during glide. Behind these obstacles the dislocations are arranged in 
active glide planes, sometimes called glide zones (Seeger 1955). The dis- 
locations in these glide zones are so far spread out (on account of the rather 
low applied stress) that they cannot be considered as piled-up groups. 
Gliding dislocations will have to pass them individually, and to this passing 
the preceding considerations regarding the effect of stress multiplication 
will be applied. 


§ 4. QUANTITATIVE FORMULATION OF THE THEORY 


Any quantitative theory of work-hardening has to proceed in three steps 
(Seeger 1956): (a) choose a dislocation model; (b) calculate the flow stress 
for this model; (c) find out how the parameters of the model vary with 
strain. Ifstep (b) has been performed, (c) gives the rate of work-hardening 
and hence by integration the stress-strain curve t=7/(e). 

The dislocation model we wish to adopt is described in §3. The experi- 
mental results suggest that to a reasonable approximation the lengths of 
the slip-lines and the density N of active dislocation sources remain 
constant. This means that the strain « and the number » of dislocations 
coming from a source should be related linearly, namely 


de _, 9 

seco Di Lig Von Te teee tee is) (ot) ee (4i) 
dn 16 a2 ; ( ) 
Here L, and Ly, are the slip-line lengths of edge and screw dislocations. 
b is the dislocation strength, and the factor (3/4)? comes in since the average 
separation between positive and negative dislocations in a slip zone is 


_ only three-quarters of its length. Expressing N in terms of the slip plane 


i 


distance x, defined in § 2.3, according to 


JERE AD 7 we maces) * Se eee ere ee SP) 
we obtain 

ei IY eas ice an A eR (4.3) 

dn 16% : 


Since x can be measured in the electron microscope, eqn. (4.3) provides the 
answer to problem (c). 

We shall now treat problem (b). The flow stress may be defined as the 
minimum shear stress at which large-scale dislocation movements are 
possible. As discussed in §3, a necessary condition for large-scale dis- 


location movement is the following: if a dislocation emitted from a 


a 


source gets stuck after a glide distance J, the back stress exerted on the 
- source must be less than or at best equal to the increase of external stress 
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during the strain interval corresponding to the emission of one dislocation 
from a source. For screw dislocations this condition can be written 


NG » Poy Cera 


The next step is the calculation of Ly. A gliding dislocation will be stopped 
if it passes too close to a parallel dislocation of the same sign. The shear 
stress necessary to push a screw dislocation past another one at the distance 
yo is equal to (Gb/27)(1/2y,). The simplest assumption we can make (which 
might not be too good if the applied stress is almost equal to the critical 
shear stress) is the following: glide can proceed if the difference between 
the stress 7 acting on the moving dislocation and the critical shear stress T» 
is at least equal to that ‘passing’ stress. 

We obtain as a second necessary condition for the continuation of glide 


Pec Reso (4.5) 
~7T “Yo 
In order to relate y, and L, to each other we have to calculate the prob- 
ability p,,, that after travelling a distance Ly a dislocation will meet 
another dislocation the glide plane of which has a distance y, from that of 
the first dislocation. This probability is determined by two events: 

1. The passing of an activated glide plane in a distance Yp. 

2. The meeting of a dislocation on this glide plane before having travelled 

a distance Rp. 

When the dislocation has travelled a distance Ro, it will meet another 
activated glide plane. This changes the internal stresses acting on the 
dislocation and the whole process is repeated. In a simplified picture we 
find (for the full calculations see Appendix B) that the probability for the 
first event is 


Py r= we 0. Me a SA SP 
and that for the second event is 
P,= — (4.7) 
2 


The combined probability p,,, is (approximately) given by the product 
of eqns. (4.6) and (4.7) 


‘ RL (4.8) 


The event we are looking for, namely the meeting of a parallel dislocation 
of opposite sign at a distance yo, has taken place if Py» lf. We obtain 
thus as the desired relation between y, and Ly + 


Rab 
Soha (4.9) 


As the condition for the flow stress we take that in eqns. (4.4) and (4.5) 
the equality signs should hold. This means that during easy glide d7/dn 


+t The possibility that the probabilities considered here may become larger 
than one is only apparent. Compare the full treatment in Appendix B. 
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and L, adjust themselves in such a way that when a dislocation meets an 
obstacle which it cannot overcome with the aid of the applied stress, the 
applied stress will be raised until another dislocation is generated at the 
source and will now help to push the previous one past the obstacle. 

By combining the preceding equations. we get (eliminating yo) 


_ 27 Rol, 7—7 (4.10) 
Uti = b G ; 
and (eliminating L,) 
bG\? n dn 
sae ee eee 
! (t7—T)) dr (=) Bola (4.11) 
Integration gives 
T=7T)+ ae G( Rei) Set est a (Ake) 
“77 


Equation (4.12) completes step (0). 
We may now obtain the expression for the stress-strain curve by com- 
bining eqns. (4.3) and (4.12). Since 


DNR, Se re fet (4513) 


we eliminate R, by using eqn. (4.2) and find as the final result for the rate 
of work-hardening 
A/G =(8/9m) . («/L,)34. SOG oe (4514) 


Equation (4.14) can be checked experimentally, since 6, x, and L, can be 
measured. On account of the great length of the slip lines, in some cases 
the experimental determination of L, is somewhat uncertain, however. 

Let us first compare eqn. (4.14) with the experimental results obtained 
on copper single crystals deformed at 90°K. In the last column of the 
table we give the values of 4; as calculated from the measured values of 
x and L,. The experimental value for the work-hardening rate of that 
crystal is 6;=0-70kg/mm?. For the other copper crystal mentioned in 
§ 2.3 (d) we calculate 0, = 0-77 kg/mm?, which is exactly the value obtained 
from the stress-strain curve. 

In Ni-20% Co alloy single crystals Kronmiiller (1959) finds x= 5004, 
L,=0-1cem. The numerical values calculated from eqn. (4.14) are at 
90°K: 0,=1-4kg/mm?, at 20°C: 6;=1-3kg/mm*. Experimentally one 
finds at 90°K: 6;=1:7kg/mm?, at 20°c: 0;=2-1kg/mm?. The agreement 
is as good as one may expect in view of the inherent difficulties of the 
problem. 

In Appendix A it is shown that the picture of work-hardening used in 
this section also explains quantitatively the results on the approach to 
ferromagnetic saturation. 

The low-temperature work-hardening of zinc single crystals requires a 
separate discussion, since ending slip lines are not or only extremely 
rarely observed. This means that the average slip distance may be 
larger than the diameter of the crystal. This situation was considered 
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elsewhere (Seeger 1958, §52) with the following result: if the number of 
obstacles for the dislocation movement is not increased during defor- 
mation (i.e. if in an infinite crystal the slip distance were independent of 
the strain), the correct formula for the work-hardening coefficient is” 
obtained by using the slip distance L,, in an infinite crystal and not the 
diameter L,, of the glide planes within the crystal. This condition is 
fulfilled for our work-hardening model. The model therefore does not 
predict a dependence of the work-hardening rate of hexagonal crystals on 
the crystal thickness, in agreement with the (rather scarce) experimental 
observations. We should get an upper limit for @ by replacing in eqn. 
(4.14) L, by L,, which under our experimental conditions was typically 
of the order 0-5em. This gives as an upper limit 06=2x10*G. The 
experimental value is indeed smaller by about a factor of two, namely 
6=10 4G. 


§ 5. Discussion 


In the present paper we give a quantitative theory for the work- 
hardening of f.c.c. metal single crystals in the easy-glide region and for 
basal glide in h.c.p. metals at low temperatures. The basic idea of the 
theory is that while the dislocations are arranged in slip zones on glide 
planes, they act through their stress fields as individual dislocations and 
not as piled-up groups. This fact is particularly strikingly borne out by 
the measurements of the differential susceptibility in the range of ferro- 
magnetic saturation in NiCo single crystals. The agreement between the 
calculated and the measured rates of work-hardening is better than one 
could have expected in view of the limited experimental accuracy and the 
statistical nature of the theory. 

The theory of stage II (Seeger e¢ al. 1957) and the present theory appear 
as limiting cases of a general theory of work-hardening, which postulates 
that the applied stress must be large enough both to drive dislocations 
through the internal stresses of the crystal and to generate new dislocations 
from Frank—Read sources against the back stresses of dislocations generated 
earlier. The present theory is the limiting case where the distance R, 
between neighbouring groups of dislocations (each of them generated from 
one source) is smaller than’ the separation between dislocations in the 
same glide plane, whereas in stage II the reverse case is realized. 

The theory developed in this paper contains two parameters that must 
be taken from experiments: the length of the slip lines (Z,) and the 
separation between slip lines (w). The situation differs therefore from that 
in stage IT, where the length and the height of the slip lines can be eliminated 
from the theory and where an expression for the work-hardening rate 
valid for all metals can be obtained (Seeger 1958). The parameters a and 
L, are determined by the dislocation structure of the undeformed crystal 
and depend therefore on the individual metal, the previous history (in 
particular the conditions of crystal growth), and the impurity content. 
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This is indeed borne out by the experiments, which show that stage I and 
the deformation of hexagonal metals are much more sensitive to impurity 
content and previous history than stage IT. 

The values of x (essentially the density of active sources) and L, (slip 
distance between the obstacles to the dislocation movement) that should 
be used in the theory may to some extent depend on external parameters. 
This has been discussed in relation to the temperature dependence of 6, 
by Seeger and Trauble (1960). The orientation dependence of 6, in f.c.c. 
single crystals (Diehl 1956a—there is no corresponding variation with 
orientation in the hexagonal metals) may have a similar explanation 
(Seeger 1957, 1958). The present theory predicts that 6; is independent 
of the crystal diameter (Seeger 1958), in agreement with the observations 
on the orientations in the middle of the stereographic triangle considered 
in this paper (Suzuki ef al. 1956). Other orientations do show a size effect 
in the easy-glide region. This can be explained by assuming that during 
deformation some new obstacles (e.g. Lomer-—Cottrell dislocations) are 
formed. A continual reduction of L, during deformation will indeed lead 
to an increased rate of work-hardening as observed for orientations 
closer to the boundaries of the orientation triangle. 

The present theory cannot hold for very small deformation. Even at 
the smallest strains investigated, a substantial number of dislocations per 
slip line (n= 10-20) is observed. This is incompatible with the notion of 
constant density and length of slip lines. We must rather expect that 
during an initial stage extending over a few per cent of shear strain the 
number of active dislocation sources increases with straint. For the case 
of zinc this is discussed by Seeger and Trauble (1960). Mader (1957) 
observed such an increase in the slip line density in copper single crystals 
deformed at room temperature, the easy glide region of which may well be 
too short to display fully the stage covered by our theory. Detailed 
- observations on that initial part of stage I are very difficult to obtain, since 
the complications discussed in § 2.2 apply a fortiore. | 

We have excluded from our discussion the work-hardening of hexagonal 
metals at higher temperatures, where recovery can no longer be neglected. 
The dislocation processes in that range appear to be rather complicated. 
The reader is referred to a recent paper (Seeger and Trauble 1960), where 
a detailed qualitative picture of these processes has been given. 
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+ If at very small strains there were slip lines with 7 less than, say, n> 5, 
we could not see them until they had grown larger and we would therefore 
_ observe an increase in slip-line density with strain. 
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APPENDIX A 
_ APPROACH TO FERROMAGNETIC SATURATION 


In this appendix we give a more detailed account of our measurements 
on the approach on ferromagnetic saturation in deformed Ni-20% Co 
single crystals. As shown in detail elsewhere (Kronmiiller 1959, Seeger 
and Kronmiiller 1960, Kronmiiller and Seeger 1961) the nature of the 
long-range stress fields due to dislocations in ferromagnetic materials can 
be studied by measurements of the differential susceptibility y=dJ/dH in 
the saturation range. These are therefore ideally suited to investigate 
whether in stage I the dislocations act as piled-up groups (as they do in 
stage II) or whether they act as individual dislocations (as is supposed 
in the theory of the present paper). 
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The general quantitative expression to start from is eqn. (6.16) of 
Kronmiiller and Seeger (1961). It was evaluated in that paper for piled- 
up groups of n dislocations. The result for n individual dislocations can 
be obtained from eqn. (6.16) by putting the distance R,, between neigh- 
bouring dislocations equal to the separation Ry between dislocation 
groups of both signs. The physical basis of this is that the stress field of 
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a dislocation is screened within the radius R,. Assuming for simplicity 
rectangular dislocation rings, we obtain 

nG?N Lb? 

4nl , 

In eqn. (Al) J, denotes the saturation magnetization, «=(HI ,(C)2 the 
reciprocal of the exchange length, and C the exchange stiffness. Equation 
(A 1) predicts about n~! times the effect of piled-up groups of n dislocations. 
With the magnitudes of n as deduced from slip-line data, this difference in 
magnitude should be easily discernible experimentally. Figure 6 gives 
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d(xH?)/dz versus 7, derived from fig. 6. The straight line is obtained 
theoretically from eqn. (A 3). Temperature of deformation: x 300°K; © 90°K. 


our experimental results for yH* as a function of the flow stress for 
H =2000 Oe and for two temperatures of deformationy. From fig. 6 we 
determined AyH?/Az, i.e. the change of susceptibility during a strain inter- 
val divided by the corresponding increase in flow stress (fig. 7). Making 


use of the relation (see eqn. (4.1)) 
16(t¢ —Te") 


nbL,N = PLS 
we obtain ae 
4 
Bier 6-6 x 10-®log («R,/2). . . (A2) 
dy [dt Oe 1 E60; x g (k« o/ ) 


Ee 
+ We plot only that portion of yH* which is due to plastic deformation 
(Kronmiiller 1959). 
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Inserting «-1=1304 and the experimental results L,=0-1cm, 
6, =2:1kg/mm?2, and Ry=1-6 x 10-*em, we finally get 
dies) 21-0852100 ee 

dtg | 1=2000 0e kg /mm* 

Comparison between experiment and theory in fig. 7 shows that the agree- 
ment is rather good. A theory based on pile-ups would have given a 
much too large effect and a wrong dependence on the applied stress. For 
this reason Kronmiiller (1959) plotted yH? against the strain ¢ in the easy 
glide region rather than the stress. Since he used (apart from a numerical 
factor) the same formula for the strain as in the present paper, he got 
satisfactory agreement with experiment. 


APPENDIX B 
CALCULATION OF SOME PROBABILITIES 


Let us first calculate the probabilities p, and p, defined in $4. It is 
convenient to divide an area A of the crystal into A/y,” ‘cells’. If Ry is 
the separation between active dislocation sources, we have to consider 
A/R,? active glide planes. As mentioned in § 4 the first prerequisite in order 
to stop a dislocation is the meeting of an activated glide plane in a distance 
Yo. There are A/&,? (the number of activated dislocation sources) possible 
events which are distributed at random over A/y,” cells. The probability p 
that a cell is occupied by such a possible event is y)?/R,?, provided 
p<i. We now wish to find out with which probability p,(Z,’) a dis- 
location gliding through a distance L,’ passes at least once within the 
distance y, from an active glide plane. This means that we have to 
consider m= L,'/yo specified cells and find the probability that not all of 
them are empty. The probability that a specific one of these cells is 
empty isg=1—p. The probability that all these m cells are empty is q™. 
The desired probability is therefore given by 


Pp, (Lo) =1—-(1—p)™. 2 5 te 


By the same kind of limiting process that leads from the binominal distri- 
bution to the Poisson distribution we obtain for p<landm>1 


Py (Lo) =[1—exp(—mp)]=mp+.... Ai 9 ea ee 


The approximation p,(Lo')=mp=L,'yo/R,2, which we will use sub- 
sequently, has a simple geometrical meaning and can easily be derived by 
inspection. 

We have now to determine the probability p, that a gliding dislocation 
actually meets one of the n dislocations in a neighbouring glide plane 
before it has passed that glide plane. The calculation is similar to the 


preceding one and gives as a simplified result under the condition eee <1 
_ 2nRy : 
p= ee ee ME Re ee 
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Equation (B 3) can also be obtained by the following simple consideration : 
the average linear density of dislocations of a given sign in a slip line of 
length L, is n/(L,/2), where n is the number of dislocation rings given out 
by the source of the line. The average linear separation between sources 
and therefore also lines (characterized by their centres, say) is Ry. The 
probability of meeting a dislocation in such a line before a new event in 
the trial leading to eqn. (B 1) takes place is therefore given by eqn. (B 3). 

Finally, we have to consider the probability Py,» that both events 
considered above take place simultaneously at least once while the 
dislocation moves a distance L,=kL,’. Py,» can again be found from 
the trial leading to eqns. (B 1) and (B 2) with 


P=Py (Lo )Py eee evi eaae (15.4) 
and m=k=L,/L,’. peas at Saheb) 
The final result is InL 

Py, o=1—exp(— a). Siravas (86) 
Cana 


We are interested in the value LZ, which makes eqn. (B 6) approximately 
equal to one. For practical purposes we may take 


Rolin = 2nL yo; Ge ME Stee ie. ber Baa es (Eig) 
which is the result used in § 4. 
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ABSTRACT 


The distribution function describing the development of electromagnetic 
cascades has been derived under the assumption that the primary particles 
are y-rays with a certain energy spectrum and angular spread. By comparing 
experimental results on the lateral distribution of electrons in cascades 
initiated by nuclear interactions with the calculations presented here one 
may obtain information on some characteristic features of ‘jets’. At large 
distances from the cascade ‘ core ’, the electron distribution is identical with 
that derived from ordinary cascade theory, thus indicating that the angular 
spread of the primary y-rays can then be neglected. This result is of interest 
when measurements in large air showers are compared with theory. 


§ 1. InTRODUCTION 
In recent years, large assemblies of particle detectors and layers of 
heavy materials have been frequently used to study nuclear interactions 
of high energy; see, for example, Proceedings of the Moscow Cosmic Ray 
Conference, 1960, Vol. I, the papers by Babayan et al., Duthie e¢ al., 
Fujimoto et al., and Grigorovy et al. By measuring the development of the 
electromagnetic cascade created by these interactions, one can obtain 
information about the energy liberated in the primary and secondary 
interactions, and of other parameters characteristic of nuclear interactions. 

If one employs counters together with dense materials, it is compara- 
tively easy to study experimentally the longitudinal development of a 
nuclear cascade. These results can then be compared with theoretical 
predictions of the development of nucleonic cascades; see, for example, 
Rozental (1952), Ueda and Ogita (1957, 1958), Fukuda et al. (1959), Olbert 
et al. (1957). The disadvantage of counters is their low resolving power 
and the complicated experimental equipment, which cannot easily be 
carried to high altitudes. 

On the other hand, the use of nuclear emulsions allows a careful study 
of the ‘core’ of the cascades; since most of the energy is carried in the 
core, it is there that most of the information is to be found. The principal 
difficulty is to relate closely the quantities that can be measured accurately 
and with confidence—such as the lateral distribution of the cascade 
particles (mainly electrons) and the variation of this distribution with the 
depth of penetration of the material in which they develop—to the 
important parameters of the primary nuclear interaction such as the 
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spread in angle and energy characterized by the transverse momentum, 
or the energy spectrum and total energy of the secondary particles. 

The aim of this paper is to take the first steps towards establishing such 
a close relationship. Since we are not interested in testing any particular 
theory of meson production, a model of the primary interaction can be 
employed which displays clearly the parameters to be determined by 
measurements on the ensuing electromagnetic cascade. 


§ 2. Tue Mops. or THE NUCLEAR INTERACTION 


Recent experimental results (Edwards et al., 1958, Co-operative 
Emulsion Group in Japan 1958) suggest that the transverse momentum 
of a particle created in a nuclear interaction is nearly independent of its 
direction of emission, and of its energy. Further, the average value, r, of 
this quantity appears to vary very slowly, or not at all, with the energy 
of the primary particle initiating the interaction (Edwards e¢ al., 1958, 
Aly et al., 1959). We therefore assume that the distribution in transverse 
momentum, p;, may be represented by the relation 


2 
Nips) drv= $25 prexp (— 2S) dpe . 3 


Figure 1 shows the distribution in p; thus determined together with the 
corresponding experimental results for pions as measured by Edwards et al. 
(1958). It may be seen that eqn. (1) is in satisfactory accord with experi- 
ment. 

We have next to consider the energy spectrum of the secondary particles. 
This may be assumed to be given by a power law 


N(E)dE cE-*dE 


in the centre-of-mass system of the interaction. Such an energy distri- 
bution is also reflected, however, in the energy spectrum of the secondary 
particles in the laboratory system, if particles projected forwards are 
alone considered. ‘Transforming to the laboratory system the energy 
spectrum of particles that are emitted forwards in the centre-of-mass 
system, and assuming that for each particle the relation »,=7 holds, one 
finds that their energy spectrum is well approximated by the same power 
law distribution with a lower cut-off at about 27y¢, were ye is the Lorentz 
factor of the centre-of-mass system. In what follows, it will therefore 
be convenient to confine attention to particles emitted in the ‘forward 
cone’ of the jet. This restriction is justified, because the particles 
emitted backwards are of very small energy in the laboratory frame of 
reference. Under these assumptions, the function describing the forward 
cone in the laboratory system of a nuclear interaction becomes 


292 
N(E, 0) dE dé = YE? Gexp ( z a *) dEd, 1<a<2, (2) 


72 


where H and @ are, respectively, the energy and tangent of the angle of 
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emission with respect to the primary direction. Equation (2) has been 
obtained by substituting HO for p; in eqn. (1), and multiplying it by the 
energy spectrum. ‘T’he normalization factor Y is obtained by considering 
the energy radiated, y: 


Eee ey pt Al en iad i sens (3) 


Number of Events 
fo} 


Pr (Be V/c) 


. The experi- 
Plot of the number of events vs. the transverse momentum p; 
ae mental results are from Edwards et al. (1958) for pions. The curve 
represents eqn. (1) with a value for 7 of 0-4 Bev/c. 


Here, the contribution from the lower limit of the integral has been 
neglected on account of the subsidiary condition of eqn. (3). 
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The transverse momentum p; may be expressed by H@ only if @, the 
tangent of the angle of emission, issmall. The lower limit of # is 27y¢> 27, 
however, so that the exponential is practically zero for values of O> 1/yc. 
It is therefore permissible to let 6 range from zero to infinity without 
introducing a great error, provided H > 2ryc. 

In the energy spectrum, the upper limit, p, has not been identified with 
the primary energy H, of the interaction, which appears only in ye. 
This has been done to illustrate the effects of such a cut-off, since doubts 
have been raised by Duthie et al. (1961), as to whether, for primaries of 
very high energy, the energy radiated in interactions caused by nucleons, 
y, is still proportional to the primary energy £). 

If «=2 (as in Heisenberg’s theory), Y (y;p) becomes 


—1 
¥ (xp) 35 (In 2) 1 2 


§ 3. DERIVATION OF THE DisTRIBUTION FUNCTION OF ELECTRONS 


In this section we derive a formula for the lateral distribution of 
electrons in a cascade which originated in a nuclear interaction of which 
the secondary particles are characterized by eqns. (2) and (3). A similar 
attempt has been made by Hasegawa (1958) who did not, however, 
give an exact distribution covering all distances from the main axis of 
the cascade. 

The procedure adopted here is to consider projected distributions 
first and then to transform to the true spatial distribution. This will 
be done in the following way. Consider the Fourier transform F'(n,, n») 
of a spatial distribution f(x,,7), where x, and 2, are the Cartesian 
coordinates in the ‘target diagram’. The relation between F and f is 


1 oe E : 
En; Ns) = = {| dx, dt, f (2,2) exp (1n,7,+inw,). . . (5) 


Let f(z,) be the projected distribution 


Flay)= |" dagflesa,), i> Saae 


==) 


and F(n,) its Fourier transform, 


Flm)= (2m) |" deyfledexp linn), . . . . (2) 
then we find the relation 
F(n,) = (27)? F'(n;, n, = 0). Soe" aan seh as Sa 


In the course of the calculation with projected distributions we will 
arrive at an expression for the function F(n,). We know that the final 
solution must be symmetrical with respect to x, and a, hence F'(n,, ng) 
must be symmetrical with respect to m, and ny. Therefore we will, 
according to eqn. (8) enlarge F(n,) in a symmetrical way so as to find a 
function F(n,, 2) which satisfies eqn. (8). This can be done in one way 
only, and it provides us with the Fourier transform of the final solution. 
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The projected distribution describing the nuclear interaction is obtained 
from eqn. (2) as 


299 
n(H, 5)dE dd = Y'(y, p)E!-*exp (-== Bide 


r") aB as, ay) 


with Y"(x, p) = (7/7) Y(x, p). 
In eqn. (9), 5 is the tangent of the projected angle of emission. 

We need in addition the projected distribution of the electron density 
in a cascade originated by a single y-ray of energy Wy. To simplify 
calculations, we consider Landau’s approximation of Kamata and 
Nishimura’s equations (1958). From their functions (3.32) together 
with (3.42), we find by using an integral representation of the Beta 
function : 


fe on ae en W\* « (ely|\-2 
TIP"! (Wo, 9, y, t) = zat. uel) a a.) 


x D(p+ 1/2)r4?1(s + 2p) 9.(p, —s— 2p, s, t) (10) 


where the notations are as in Kamata and Nishimura’s paper, and y is 
the projected distance from the cascade core. For the purpose of 
simplifying the notation we introduce 


ee =o ss) Wyler e ds. eee) 
with the integral operator 


l +40 —2p B : 
$(p,s)= — zal{_,. ds dp —(F) D(p + 1/2) 0(s + 2p) at? 


x NUP, —s—2p,8,t). oR es (12) 
The distribution of electrons is found as follows. In the primary 
interaction, a number 


n, AE d= Y'(y,p)E'* exp(- 


pee 


*) abs pera (13) 


of neutral pions are emitted into dH dé at H and 6. In eqn. (13), and 
in what follows, the quantity y describes the energy radiated into the 
neutral pion component. Hach 7°-meson decays into two y-rays almost 
immediately. The probability distribution of finding a y-ray at energy 
W, in dW, is well approximated by a rectangular distribution, extending 
from 0 to #, 

VET ag <li! 2 Gen 7 oes eee (14) 


The maximum transverse momentum of a y-ray with respect to the 
direction of its parent 7°-meson is 70Mev/c, much smaller than 7. Thus 
we can neglect the angular divergence of the y-rays from the direction 
of motion of their parent pion, and assume them to be collinear with it. 
Hence the number of y-rays emitted from the jet into dW, do at Wy 


and 6 is 


p 1 
n,dW,d3=2dW ,d5 | dE n,., (15) 
1) 


E’ 


_where the lower limit 6 has the value W, or 2ryc, whichever is the larger 
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Each of the y-rays creates a cascade. Consider the electron density 
at a projected distance a from the core, and a distance ¢ from the jet 
(all distances are measured in units of the radiation length). According 
to eqn. (11), each y-ray of energy W, contributes a density 


d(p,8)W y*|st—a|-"P9 de. se 2 ee KEG 


towards the total electron density D(x) at x. This total electron density 
is found by integrating over all angles of emission, 6, and all photon 
energies, Wo: 


D(a) =d3 |" dWoplp,8)W.2[8t—2h In, ie ca EGE 
0 


The integral over dS can be solved by a Fourier convolution (Sneddon 
1951). 


+ 00 2 
| d8|St — a|-2P 4 exp(— mbes! 


)- (ioe dn exp [—in(x/t) ] 


x A(n) Blnjt-? =. es 
with 


J ity\8 |) dig lal eet ot te T'(— 2p)nir 
A =| 27T [2 p = qT 
(n) = (27) [oa |&|-??-1 exp (ink) = (27)? Tpe+pFd2op)’ 


he ‘ E?5°x 2\12 n*7* 
dk ex k— ay og Yee 
a exp (in nz ) (=) Bi exp (-5 Ba 7): (19) 


One obtains the following expression for the electron density: 


B(n) = (20) | 


dx 


D(wyde= {: dW b(p, 8) Wer P(x, p |" dE R-= | fea exp| —in() | 


nr D'(— 2p) 
x exp( — <> )\t-% 0? aie 
(-3 Ber ) D(1/2+p)P(1/2—p) a 
At this stage, the transformation to the spatial distribution can be 
performed. Identifying n with n, and x with x, in eqn. (20), one realizes 
that the Fourier transform of the projected distribution, F(n,) is 
2-2 


“ > ‘ 
F(n,) = (2m)! | set (p, 8) W847) (x) |" dE E---1 exp(- ie ) 


r(=2p) 
x (277,27 Ninn ee 2 

TOR + pr a2=p) eT 
Therefore, the double Fourier transform of the true spatial distribution, 
F(ny,N) must be, according to eqn. (8), 


F'(n4, Ng) = |; dW h(p, 8)W 47 Y'(y, P| . dE B-exp ( e ae 
b 7 


{-2P(m,.2 2 I'(— 2p) 
x ENP(ny +")? Tapp) Rop)’ sy (2a 


since it yields eqn. (21) for n»=0 and it is symmetrical in MN, and np. 
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The spatial electron density distribution is therefore 
Ye, ts)d0 dt. 


dx, dx + ae ; 
= a = (2Qzr)-1 [| dn, dn, Fn, 2) exp (s ny 2 —iNg 2) : 
“ =O) A / 


(23) 
Introducing r= (x,?+2,7)1?, m=(n2+n,2)"? and v, the angle between 
the vectors (%,%) and (n,,.), we obtain by transforming to polar 
coordinates : 


oD 


p 
Dr) =| dW $(p,9)Wo'4r¥'(x,p)| ABE (ny dmm m?? 
0 b 


0 


Qn — 22 T'(—2p) 
dvoxpl 2am NeT ig \ Ey cae ey a ct 2) 
x iE vexp( um. ; COS v a Pdj2+prapos): (24) 


This turns into the inverse Hankel transform: 


D(r)= [ arrose. a) Worse ¥"(x,) | amu)" dmm me. J,(m- ) 
0 , b 5 
Bea) pee yeareee Sk AP) a 
x exp ( i) T(ij2+p)r(j2—p) i (25) 


For the solution of these integrals we introduce the Mellin transform 
with respect to H, i.e. we write: 


m7? 


+40 
’ IE B-=texp(~ )- | ° dE(2ni) | duE-"1M(u), (26) 
b E TT b —4o 


with 
et heey aN mt \¥* 1, fo —u 9 
M(n)= | dE E exp ( Te) =12(Sr) r( 5 s (27) 


In doing this, the integrals over dW, and dH and the inverse Hankel 
transform can be solved. The final solution is: 


+400 ad as” 
D(r) or as: (isnt | ff ds dp du ¥'(x.p)2e( 3] (Fr) 


s 


a—2p—u—2 
x (5) ‘ tp? 2D (s + 2p) 1 — 2p) 2(p, Rie 8, t) 


2t 
p (Pee) 0") é 
eee pee eee eS SE (pe Ot es (2 s—utl)_ 
Be (4) *) aengcaey. (ry 0) a 
2 ota 


Fee Cee) 


§ 4, DISCUSSION OF THE SOLUTION 


If one wants to consider the distribution of electrons for small values 
of r, one may follow the method employed by Kamata and Nishimura 
(1958) in a similar situation. This is to consider, in the integral over p, 
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the contribution from the pole at p=—s/2 of I(s+2p) only. Then, 
eqn. (28) becomes 


+400 u—% py \ats—u-2 
D(r)= — (isn) | . ds du Y'(x, p)27 (=) B=) 


1 
(st+t1)(st+1l—u 


(eS) FS) 
ee eS Mok) pac let Eon ee oe ae 
GEE) 
2 2 


For very small r, the pole at ~+2—a—s=0 is leading, so that the 
distribution is flat in the middle. Furthermore, s+1—u=3—a>1, so 
that the term with p is predominant. It should, perhaps, be emphasized 
that an electromagnetic cascade initiated by a single particle has also a 
flat distribution in the middle, just like the function of eqn. (29). This 
is due to the fact that one can write the electron density for small values 
of r in a cascade produced by a single particle in the form 


Il, nes) ps—2 


x -ntl2D(s) $92 — 4/2, 0, 8,0 (oe tS Baye) 


where s is determined by the equation 


Ai (s)t-+In 2” —1)24(=5*) =0. eS eae 
E, 2 

Equation (30) is slightly different from the corresponding equation 
(3.46) of Kamata and Nishimura (1958) where the % function has been 
dropped. Though this is permissible for larger r, this is not the case 
when r is small, since it is the % function which compensates the large 
negative values of the logarithmic term, so that s+2 for r+0 according 
to eqn. (30). Therefore 


Il,=const., s=2 


in the centre of a single electromagnetic cascade. 
The integral number 
R 
N(R) =2n | rDir\de 

0 
of electrons within a circle of radius R has been calculated by a double 
saddle-point method, and numerical values have been obtained for the 
special assumptions 


p=, and l0y=F). 


It has to be remembered that these results have been obtained under 
Landau’s approximation of the cascade function. However, the results 
can be normalized to the cascade function of a single photon in Landau’s 
approximation. This normalized distribution can finally be multiplied 
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WoR ond xR respectively (units BeVu,rad.L= 2.9cm) 


Plot of the integral number of electrons N(f) vs. the quantities W, times R 
and y times # respectively. Curve A is the distribution of electrons 
caused by a y-ray incident at zero radiation lengths, according to 
Nishimura (1960). Curve B is the distribution due to a jet, when an 
energy y was radiated in form of 7°-mesons with an average transverse 
momentum 7=0-4 Bev/c, and an energy spectrum with an exponent 
a=1:5. The inelasticity y/H, was assumed to be 1/10. The graph is 
for pure emulsion with a radiation length of 2:9cm. The distance to the 
jet or origin of the y-ray is ¢=8 radiation lengths. 
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by the true cascade distribution obtained without using Landau’s 
approximation, and the result is displayed in fig. 2. Curve A represents 
the number of electrons vs. the quantity W Rk (W, is the energy of the 
primary photon), curve B is the distribution due to a jet with 10y=£p, 
and the bottom scale is in units of yf. 
The meaning of the term 
1 s—ut+1 9 s—u+l1 
DU tttttEEtIEE EES U — ({2T' 

in eqn. (28) can be elucidated in the following way. Consider a pencil of 
y-rays incident upon an absorber. Their spectrum may be 


E-“dE, 


and it may have the upper and lower limits p and 2ry, respectively. The 
electromagnetic cascade created by these y-rays can be described by 
the well-known formulae for the cascade development (Rossi 1956, 
Kamata and Nishimura 1958), where the term W,* containing the primary 
photon energy has been replaced by A(s,u). This quantity therefore 
represents the energy term of the ‘collinear approximation’ of all the 
y-rays from a jet. Since w varies with r, however, the effective slope of 
the energy spectrum, w, is not constant with the radius. 

For larger values of 7, the pole at ~a=w in eqn. (28) becomes leading. 
Evaluating eqn. (28) at this pole, one finds: 


Dir)==— (dney[ [ ds dp27 Y'(x, p)B(s, «)e8(e/H 221-222 


x T(s+2p)T(p+ I) galp, —8—2p,8,t), . . . . (31) 


This is identical with Kamata and Nishimura’s (1958) expression (3.28) 
for the lateral distribution of electrons, apart from the energy term 
which is now A(s,«). This means that, away from the core, the shower 
becomes identical with one being created by a narrow collinear beam 
of y-rays, with energy spectrum and cut-offs according to the jet model. 
Therefore, one may use the distribution functions of Kamata and 
Nishimura (1958) for larger values of r, ie. their expressions of §3.C, 


together with their figs. 2 to 5. The shower age s is determined in a 
slightly different way, like 


Atta) din Bie ye ee (32). 
Os 
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ABSTRACT 


The lattice spacings of h.c.p. ¢ phases in the systems Ag-Cd, Ag-Hg, 
Ag-Ga, Ag—In, Ag-Sn, Ag—As and Ag-Sb have been studied by the x-ray 
powder method. 

The results show consistent trends in the changes of lattice spacings when 
plotted as functions of electron : atom ratio. These trends may be inter- 
preted in terms of distortions of the Brillouin zone following interactions 
between the Fermi surface with various zone faces. The results indicate 
that in Ag-based alloys the overlap of Fermi electrons across the {10-0} faces 
of the Brillouin zone occurs in the region of e/a between 1-39 and 1-40. , 

The influence of particular solute elements on the changes in lattice 
spacings may be interpreted qualitatively by considering that the band 
gaps in the Brillouin zone of silver are decreased by the addition of the 
solute element. 


§ 1. INTRODUCTION 


INTEREST in the theory of the properties of h.c.p. metals and alloys was 
initiated by Jones (1934) when he applied the Bloch Theory to interpret 
the lattice spacings of the « and 7 phases in the Cu—Zn system. However, 
it was not till some twenty years later that his theory was used to discuss 
the occurrence and stability of the h.c.p. ¢ phases in other noble metal 
alloys (Raynor and Massalski 1955, Massalski 1956, 1957, 1958), and 
the lattice spacings of the latter phases were studied in detail. The 
¢ phases, which occur in alloys of the noble metals with the B-sub-group 
elements of the Periodic Table have been included in the group of inter- 
mediate phases classified as ‘electron compounds’ (Hume-Rothery and 
Raynor 1954). These phases usually follow the primary solid solutions 
in the phase diagrams of the respective alloy systems and often extend 
over wide ranges of electron concentration{ within the general range 
between 1-2-1°8. 

In recent studies the ¢ phases in two systems based on copper 
(Massalski and Cockayne 1959) and in four systems based on gold 
(Massalski and King 1960a) have been examined in detail. Several 
consistent trends in the changes in lattice spacings and axial ratios were 
recognized when these were plotted against the electron concentration. 


+ Communicated by the Authors. 
{ The term electron concentration, or ¢/a, is used to denote the ratio of all 


valency electrons to the number of atoms. 
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Additional interesting information was obtained by comparing the lattice 
spacings of the ¢ phases with those of the respective primary solid solutions. 
The tentative conclusion reached was that the lattice spacings are influenced 
by interactions between the Fermi surface and the Brillouin zone. The 
trends in lattice spacings observed for the ¢ phases based on gold differed 
significantly from those for the copper-based. alloys. Therefore, it was 
proposed (Massalski and King 1960a) that the Fermi surface at the 
beginning of the gold-based ¢ phases is probably considerably distorted 
from the spherical shape; whereas in alloys based on copper it is likely 
that the Fermi surface is nearly spherical. Furthermore, the degree of 
distortion of the Fermi surface was considered to be somewhat greater 
in alloys of gold with Hg or In than in gold-based alloys containing Cd 
or Sn. These conclusions were shown to be in general agreement with 
the model suggested by Cohen and Heine (1958) to allow for the influence 
of the solute element on the band structure of a noble metal solvent. 

Following the studies of copper and gold-based alloys, and the con- 
clusions reached regarding their band structure, it thus became of interest 
to investigate in a similar manner the lattice spacings of the seven ¢ 
phases formed when Ag>is alloyed with the following B-sub-group 
elements of the Periodic Table: 


— Ag-—Ga — Ag—As 
Ag—Cd Ag—In Ag-Sn Ag-Sb 
Ag-Hg — — — 


The materials and methods used in this investigation are described in 
the Appendix, § 1. 
§ 2. RESULTS 


The values of the a and c lattice spacings and of axial ratio, together 
with the composition and other pertinent data, of the alloys studied are 
listed in the table. In accordance with a preliminary metallographic 
examination, diffraction lines of a second phase, in addition to the h.e.p. 
pattern, were observed in alloys Ag—Cd 26, Ag—Hg 1, Ag—Hg 6, Ag-Ga 6, 
Ag-In 3, Ag-Sn 13 and Ag—As 1. The electron: atom ratio values for 
these alloys listed in the table were estimated by locating the vol/atom 
data on the respective plots in fig. 1. A comparison between the present 
results and those reported previously in the literature is given in the 
Appendix, § 2. 

The lattice spacings and volumes of the unit cell per atom (vol/atom) 
for the silver-based alloys are plotted in fig. 1 as functions of electron 
concentration, e/a. This figure has been arranged to follow the plan of 
the Periodic Table for the B-sub-group elements which have been alloyed 
with silver. In its own second long period silver forms £ phases with all 
the elements from Groups I[B-VB, and with the exception of Ag—Cd 
these phases extend over wide ranges of electron concentration, The 
absence of ¢ phases when Ag is alloyed with Tl, Pb or Bi from the third 
long period may be due to the large differences in atomic size between 


Composition, heat-treatment, and lattice spacings at 30°c of h.c.p. 


¢ phases in silver-based alloys 
eee ee ee ee eek ORS ee 


2°9276 
2:9327 
2-9368 
2-9439 
2°9553 
2-9616 
2-9646 
2-9681 


4-7831 
4-7894 
4-7920 
4-7950 
4°7973 
4.7990 
4-7997 
4-8003 


+ e/a estimated by locating the vol/atom on the respective plots given in fig. 1. 


Heat- 
treatment l 
At. % | At. % ‘ Vol/ 
Alloy ees we Structure Cn (Om UN c/a atom, e/a 
Temp. | Time, As 
eG days 
Ag—Cd 
23 49-28 | 50-72 438 60 C 3-9867 | 4:8253 | 1-6156 | 18-638 | 1-507 
24 48-18 51-82 ra 2:9904 | 4-8244 | 1-6133 | 18-681 | 1-518 
25 46-11 53-89 t 2°9976 | 4-8237 | 1-6092 | 18-768 | 1-540 
26 44-87 | 55-13 C+y 2-9983 | 4:8234 | 1-6087 | 18-776 | 1-542 
Ag—Hg 


1 | 60-00 | 40-00 | 150 | 100 a+ | 2-9758 | 4-8417 | 1-6270 | 18-565 | 1-418+ 
3 | 56-50 | 43-50 t 2-9853 | 4-8413 | 1-6217 | 18-683 | 1-435 
4 | 55-50 | 44-50 t 2-9896 | 4-8413 | 1-6194 | 18-736 | 1-445 
5 | 54-00 | 46-00 : 2-9990 | 4:8410 | 1-6143 | 18-853 | 1-460 
6 9.5 aon €+second Vins ees : : 

52-50 | 47-50 phase | 20004 | 4-8408 | 1-6134 | 18-870 | 1-463+ 

Ag-Ga 
2 |78-30|21-70| 462 44 t 2-8820 | 4-6933 | 1-6285 | 16-880 | 1-434 
3 | 77-50 | 22-50} 540 21 t 2-8827 | 4-6896 | 1-6268 | 16-875 | 1-450 
4 | 76-50 | 23-50} 600 14 2-8840 | 4-6851 | 1-6245 | 16-874 | 1-470 
5 | 75-00] 25-00} 575 14 t 2-8869 | 4:6753 | 1-6195 | 16-872 | 1-500 
6 | 73-50] 26-50] 540 a1 | §+second | 9. g295 | 4.6645 | 1-6143 | 16-864 | 1-528t 

phase 

Ag-In 
3 | 78-00 | 22-00] 540 21 att | 2-9455 | 4-7912 | 1-6266 | 18-000 | 1-448} 
4 | 76-50 | 23-50} 655 14 " 2-9495 | 4-7891 | 1-6237 | 18-041 | 1-470 
' 5 | 75-00 | 25-00 | 600 14 c 2-9563 | 4-7857 | 1-6188 | 18-111 | 1-500 
6 | 73-50 | 26-50} 600 14 E 2-9642 | 4-7821 | 1-6133 | 18-194 | 1-530 
7 | 72-00 | 28-00] 600 14 t 2-9739 | 4-7767 | 1-6062 | 18-293 | 1-560 
8 | 70-50} 29-50] 540 21 t 2:9840 | 4-7690 | 1-5982 | 18-388 | 1-590 
9 | 69-00] 31-00] 540 21 2-9950 | 4-7585 | 1-5888 | 18-483 | 1-620 
10 | 77-50 | 32-50] 462 44 t 3-0054 | 4-7488 | 1-5801 | 18-573 | 1-650 

Ag-Sn 
2 | 87-10]12-90]| 395 21 t 2-9284 | 4-7844 | 1-6337 | 17-766 | 1-387 
3 | 86-67 | 13-33 2-9302 | 4-7850 | 1-6330 | 17-790 | 1-400 
4 |86-00|14-00] 475 t 2-9335 | 4-7854 | 1-6312 | 17-832 | 1-420 
5 | 85-00 | 15-00 t 2-9386 | 4-7860 | 1-6287 | 17-896 | 1-450 
6 | 84-00 | 16-00 £ 2-9436 | 4-7845 | 1-6254 | 17-951 | 1-480 
7 | 83-00 | 17-00 £ 2-9491 | 4-7840 | 1-6222 | 18-017 | 1-510 
8 | 82-00} 18-00 2-9548 | 4-7835 | 1-6189 | 18-084 | 1-540 
9 | 81-00 | 19-00 i 2-9598 | 4-7822 | 1-6157 | 18-141 | 1-570 
10 | 79-99 | 20-01 t 2-9658 | 4-7824 | 1-6125 | 18-215 | 1-603 
11 | 79-60 | 21-40 t 2-9730 | 4:7818 | 1-6084 | 18-301 | 1-642 
12 | 78-00 | 22-00 2-9760 | 4-7824 | 1-6070 | 18-342 | 1-660 
13 | 77-50 | 22-50 | - ¢+second | 9 9783 | 4.7819 | 1-6056 | 18-367 | 1-672+ 

phase 

Ag—As 
1 |90-50| 9-50| 540 21 a+¢ | 2-8985 | 4:7347 | 1-6335 | 17-224 | 1-400f 
2 | 89-50 | 10-50 t 2-8993 | 4-7346 | 1-6330 | 17-233 | 1-420 
3 | 88-50 | 11-50 t 2-9008 | 4:7347 | 1-6322 | 17-252 | 1-460 
4 | 87-50 | 12-50 t 2-9025 | 4:7339 | 1-6310 | 17-269 | 1-500 
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the solvent and solute in these systems, a measure of which He bai 
by the size factors (Hume-Rothery, Reynolds and Raynor i she 
at the bottom of fig. 1. This is consistent with the bead tet a is 
ranges of homogeneity of the primary solid solutions of Ag wi ; 
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Pb and Bi are more restricted (fig. 1) 
Ag-based systems. 


Electron Concentration, e/a 


Variation of lattice spacings with electron concentration 
in-a and ¢ alloys of silver with the B-sub-group elements. 


compared to those for the remaining 
The size factors for the systems Ag—Ge 
are on the borderline for the formation of h.c 
Reynolds and Raynor 1940). 

phases under equilibrium c¢ 


and Ag—As 


-p. € phases (Hume-Rothery, 
Attempts to discover intermediate electron - 


onditions in Ag-Ge have so far proved 
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unsuccessful (Hume-Rothery, Raynor, Reynolds and Packer 1940). 
However, Duwez ef al. (1960) have recently reported that metastable 
h.c.p. phases may be induced in this system by very rapid quenching 
from the melt. It is also possible that electrochemical effects (Hume- 
Rothery, Reynolds and Raynor 1940) play a role in limiting the range 
of homogeneity of the ¢ phase in the Ag—As system and possibly also 
in Ag-Hg. The tendency towards ordering at low temperatures stabilizes 
the b.c.c. 8 phases in the systems Ag-Cd and Ag-Zn at the expense of 
the disordered h.c.p. ¢ phases, even though the size factors appear to 
be favourable for the latter phases. In Ag—Cd the onset of the ¢ phase 
is displaced towards higher contents of cadmium, and its range of 
homogeneity is further limited by the presence of the more stable y 
phase in the region of e/a from 1:55-1:63. In Ag—Zn the tendency towards 
order is such that a pseudo hexagonal phase, £°, is formed in place of the 
usual ¢ phase, but a h.c.p. structure may be induced by cold-working 
the ordered f’ phase (obtained by quenching f) (King and Massalski 1961) 
at room temperature. It is of interest to note that at values of e/a 
beyond the y-phase range (i.e. for e/a > 1-7) extensive h.c.p. € phases are 
observed in both Ag—Zn and Ag—Cd. 

As in the previous study of gold alloys (Massalski and King 1960a) 
values of a’ and hypothetical c’ spacings and of vol/atom for the « 
primary solid solutions have been included in fig. 1 for comparison with 
the present results. The values for the «-phase alloys were calculated+ 
from the data of Owen and Roberts (1939) for Ag—Cd, Ag—In, Ag—Sn 
and Ag—Sb; the data of Owen and Rowlands (1940) for Ag—Ga and Ag—As; 
and that of Day and Mathewson (1938) and Rayson and Calvert (1958) 
for Ag—Hg. 

On passing from the primary solid solution to the ¢ phases in the 
systems Ag—-Sn, Ag-As and Ag-—Sb there is a small but measurable 
decrease, Aa’ (fig. 1) in the a lattice spacings, and a complementary 
increase, Ac’, in the c spacings. These deviations are similar in form, 
but very much smaller in magnitude, to those observed previously in 
gold-based alloys (Massalski and King 1960a). As a result of such 
deviations the axial ratios at the onset of a number of ¢ phases lie above 
the ideal value for the close-packing of hard spheres, as shown in fig. 2. 

Immediately following the onset of the ¢ phases in Ag-Sn, Ag—As 
and Ag-Sb the curves of the a lattice spacings appear to lie parallel to 
the a’ spacings curves for the respective primary solid solutions. How- 
ever, at values of e/a in the region 1:39-1:40 the a spacings show an 
expansion Aa which is increased with increasing electron concentration 
as shown in fig. 1. Although the ¢ phases in the remaining silver-based 
systems are not stable until values of e/a greater than 1-4, they also show 
Aa deviations in the a spacings with respect to the a’ spacings extrapolated 


30°c using the thermal coefficients of expansion reported by Owen and Roberts 
(1939) and interpolating for the systems Ag-Hg, Ag—As and Ag—Ga. 
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from the respective primary solid solutions. In the Ag-Ga system the 
Aa deviation is such as to cause the a spacings to increase ed 
e/a whereas the addition of Ga to silver causes the lattice spacings ss hi 
primary solid solution to be diminished. This behaviour is similar to 
that observed earlier in the system Ag—Al (Massalski and Cockayne 


1959). 


or tae 


Ideal Axial Ratio 


~ 


163 

[ Ag-Cdo. Ag-Snm 
fF AGHI® Age aso 
Lolweeatcas Sa? * 
[ 


Ds Ag-Ina Ag-Al x 


1.61 ees 


Au-Cd 0 
Au-Hg e 
Au-In & 
Au-Sn @ 


= es eal 


Axial Ratio, c/a 


o 

fo) 

ail 

o 

By 
T 


1.633 


rape fo 
as 
fo) 
WwW 
[tt 


159 


Cu-Ga A 
[ousec 


153-160 —— ee ee eee 
es 14 ES 


1.0 14 1.2 (PS) 1.4 iS 1.6 ura 18 
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Axial ratio changes with electron concentration in h.e.p. € phases based on 
silver compared with data for alloys based on copper and on gold. 


The c lattice spacings also show a decrease in slope with respect to 
the hypothetical c’ curves in the region of e/a near 1:39-1:40. The form 
of this deviation is also dependent on the trend in lattice spacings in the 
« primary solid solutions. In the Ag-Ga system for which the latter 
trend is a decrease in a’, and hence also a decrease in c’, the c lattice 
spacings decrease at an accelerated rate. In Ag—Sb and Ag—As, however, 
the rate of increase of the c spacings is merely diminished, whereas in 
Ag-Sn the ¢ spacings pass through a shallow maximum, The deviation 
is more marked in the remaining systems since, even though the c’ 
spacings increase rapidly, the trend is for the c spacings to be decreased. 
In general the degree of the change in slope of the ¢ spacings follows the 
magnitude of the Aa deviations for the respective systems, both effects 
being greatest in Ag—In and least in Ag—As. 
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The accelerated increase in the a spacings and the simultaneous decrease 
in the ¢ spacings at e/a=1-39-1-40 causes a rapid fall in axial ratio to 
values below the ideal. The axial ratios of ¢ phases in alloy systems 
based on copper (Massalski and Cockayne 1959) and on gold (Massalski 
and King 1960a) also show this rapid fall in axial ratio (fig. 2), which 
is one of the distinctive characteristics of the ¢ phases (Massalski 1956, 
1957, 1958). However, the present results for the silver-based ¢ phases 
reveal that the rate of the decrease in axial ratio with increasing electron 
concentration is not constant but varies markedly with the particular 
solute alloyed with silver, as shown in fig. 2. Nevertheless, the axial 
ratios in all systems, with the possible exception of Ag—Al (Massalski and 
Cockayne 1959), may be extrapolated to a common value ~ 1-6337 
(i.e. slightly greater than the ideal) in the region of e/a between 1-39 and 
1-40. 

The vol/atom at the onset of the silver-based ¢ phases appears to fall 
almost exactly on the respective curves extrapolated from the « primary 
solid solutions. With increasing e/a, however, the volume shows a 
departure towards higher values, and the vol/atom curves rise above the 
curves extrapolated from the primary solid solutions. At the onset of the 
¢ phases in the gold-based systems the vol/atom curves were also found 
to superimpose upon those extrapolated from the «a phase. 


§ 3. Discussion 


The significance of the lattice spacings results will be considered in 
general in terms of their changes with electron concentration (e/a). For 
this purpose we shall assume that both the f.c.c. and h.c.p. phases are not 
defect structures and that all lattice sites are occupied in each unit cell. 
Similarly we shall assume that the conventional valencies usually 
associated with the noble metals and the B-sub-group elements are 
applicable; i.e. Cu, Ag, Au=1; Zn, Cd, Hg=2, etc. 

Since the Brillouin zone of the h.c.p. structure does not possess cubic 
symmetry, interactions between the Fermi surface and various zone faces 
may result in distortions of the zone without altering its basic symmetry. 
Such zone distortions should in turn cause changes in lattice spacings 
along the related directions in the real lattice. Conversely, deviations 
observed in the trends in lattice spacings as functions of electron con- 
centration may be used to provide information concerning possible 
changes both in the shape of the Fermi surface and in the nature of the 
band structure, which occur on alloying. 

A section through the Brillouin zone proposed by Jones (1934) is 
shown in fig. 3 for an h.c.p. alloy of ideal axial ratio. This zone is bounded 
by 20 faces: six of the {10-0} type, two of the {00-2} type, and twelve of 
the {10-1} type. According to the original model of Jones (1934) overlap 
of electrons across a zone face may cause that face to be moved closer 
to the origin of k space and hence may result in an expansion of the 
corresponding lattice spacings in real space. The accelerated expansions 
in the a spacings of the h.c.p. alloys, marked Aq in fig. 1, are most likely 
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associated with such a mechanism; i.e. in this case the overlap of Fermi 
electrons across the {10-0} faces of the Brillouin zone which are nearest 
to the origin (fig. 3). The simultaneous contraction of the c spacings is 
very likely a compensating effect for the Aa expansion, and hence the 
detailed form of the c spacings curve (e.g. slope, presence of maximum, 
etc.) is of no special significance as discussed previously for Au-based 
alloys (Massalski and King 1960a). However, this compensation appears 
to be only partial, since an increase in the vol/atom is observed (fig. 1) 
as the electron concentration is increased beyond about 1-4. 

A measure of the relative expansion of the a spacings for the various 
silver-based ¢ phases may be obtained from the plots of the fractional 
distortions Aa/(a— Aa) plotted as functions ofe/a in fig. 4. On the basis 
of similar plots for the systems Cu-Ge, Cu-Ga and Ag—Al it was suggested 
that the rates of distortion were nearly identical for these three systems 
(Massalski and Cockayne 1959). However, the present results reveal 
that any such agreement between the rates of distortion in the copper-based 
alloys and in Ag—Al must be regarded as coincidental, since the rates of 
distortion vary over wide limits depending on which particular solute 
element is alloyed with Ag. Nevertheless, all the plots Aa/(a—Aa) for 
the B-sub-group Ag-based alloys extrapolate to zero at values of 
electron concentration which lie in a narrow region of e/a between 
1:39-1:40. If the above interpretation of the Aa expansion is accepted, 
these values of e/a may be taken as the points at which overlap occurs 
across the {10-0} faces in the various alloys. The equivalent ranges of 
e/a associated with overlap aeross the {10-0} zone faces in ¢ phases of the 
other noble metals were reported to be 1-415—1-:43 in Cu-based alloys 
(Massalski and Cockayne 1959) and 1:35-1:36 in Au-based alloys 
(Massalski and King 1960a). Hence, the value of electron concentration - 
associated with overlap varies considerably from one noble metal solvent 
to another, but only slightly when various solutes are alloyed with any 
one of the noble metals. An exception is the value of 1:47 extrapolated 
for Ag—Al (Massalski and Cockayne 1959) indicating that the A-sub-group 
element Al when alloyed with silver has an effect which is different from 
that of the B-sub-group elements. 

Studies of the anomalous skin effect (Pippard 1957), the de Haas-van 
Alphen effect (Shoenberg 1960) and of ultra-sonic attenuation (Morse 
et al. 1960) indicate that the Fermi surface in all three noble metals, Cu, 
Ag and Au, may be considered to a first approximation as a sphere with 
small cylinders (referred to as ‘necks’, Shoenberg 1960) which make 
contact with the {111} faces of the f.c.c. Brillouin zone. The cross- 
sectional area of these necks is considered to be greatest in copper and 
least in silver. For the case of a hypothetical ideally close-packed h.e.p. 
noble metal it is thus reasonable to assume that the Fermi surface would 
contain analogous cylinders making contact with the {00-2} faces of the 
zone, since these faces lie the same distance from the origin as the {111} 
faces of the f.c.c. zone and are related to similarly close-packed planes 


Fig. 3 


Vertical section through the Brillouin zone of h.c.p. structure 
of ideal axial ratio (after Jones 1934). 


Fig. 4 


- Electron Concentration, e/a 


Plot of the distortion due to overlap [ Aa/(a—Aa)] against electron 
concentration for h.c.p. phases based on silver. 


’ 
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in the real lattice. Considerable contact would probably also take place 
between the Fermi surface and the {10-0} faces which are nearest the 
origin (see fig. 3), but in directions perpendicular to the {10-1} faces, 
which are furthest away, it is likely that the Fermi surface would remain 
approximately spherical as in directions perpendicular to the {200} 
faces of the f.c.c. zone. 

If it is considered that the degree of contact, i.e. the size of the necks, 
is indicative of the band gaps across the zone faces, and that these band 
gaps are modified by alloying as suggested by Cohen and Heine (1958), 
changes in the Fermi surface due to alloying may be predicted in a 
qualitative manner. In Cu, for which Hp<Hs, the band gaps are con- 
sidered to be decreased by alloying. Hence, the diameter of the necks 
should be decreased while at the same time the volume enclosed by 
the main body of the Fermi surface should become increased due to the 
change in e/a. Conversely, in Au, for which it is considered that Lp > Es, 
the band gaps should increase on alloying, and the size of the necks should 
become even greater, at the expense of the volume of the main Fermi 
surface. The main body of the Fermi surface, being approximately 
spherical, will make direct contact (ie. without the need of additional 
cylinders) with the {10-0} faces of the h.c.p. zone at e/a near 1-14. How- 
ever, in order to enable the contours of constant energy to intersect the 
zone face at right angles, the Fermi surface will tend to spread along 
the zone faces resulting in a short, wide neck. An increase in electron 
concentration due to further alloying will in general cause the Fermi 
surface to spread further, and hence to a first approximation it may be 
considered that the area of contact will be independent of any changes in 
the band gaps across the {10-0} faces. Any application of the Cohen 
and Heine (1958) model to the h.c.p. € phases should thus be limited to 
the {00-2} and perhaps the {10-1} zone faces. Interpreting the lattice 
spacings results in terms of these predictions, the much higher-than-ideal 
values of e/a observed in Au-based alloys (fig. 2) appear to be in some 
way associated with a large degree of contact between the Fermi surface 
and the {00-2} faces of the h.c.p. zone. For Cu-based alloys the axial 
ratios are only slightly greater than the ideal and this, in comparison, 
would be consistent with a relatively smaller amount of contact between 
the Fermi surface and the {00-2} faces. 

If the magnitude of the departure of the axial ratio above the ideal 
value may be taken as a qualitative measure of the degree of distortion 
of the Fermi surface from a spherical shape, the Fermi surface in Ag-based 
alloys may be considered to be almost spherical in the <00-2) directions 
since the axial ratio is only just greater than the ideal, being 1-6335-1-6338 
(fig. 2). In this respect it is of interest to note that the Ag-based ¢ phases 
of almost ideal axial ratio occur in the region of electron concentration 
near 1:36, which is the value at which a spherical Fermi surface would 
just touch the {00-2} faces of the Brillouin zone, assuming any distortion 
due to interaction with the {10-0} faces to be quite local. Hence, it may 
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be inferred that the necks, which form part of the Fermi surface in pure 
Ag, are reduced in size by alloying until they tend to almost disappear. 
This indicates that, according to the model of Cohen and Heine (1958), 
Hy < Hs in pure Ag, i.e. that the band structure of silver is similar to that 
of Cu. 

Since it appears probable that in h.c.p. alloys of the noble metals the 
_ Fermi surface would remain in contact with the {00-2} faces over the range 
of e/a from 1-0 up to the onset of these phases at e/a=1-2-1-36, earlier 
interpretations (Massalski and King 1960a) of the Ac’ deviation, and 
the higher-than-ideal axial ratios, in terms of an electron stress effect 
(Goodenough 1953) as the Fermi surface approaches the {00-2} faces, must 
be reconsidered. It might be supposed that the electrons in the necks 
of the Fermi surface adjacent to the intersection with the {00-2} exert an 
attractive force on the zone boundary, and that the total electron energy 
may be lowered if the boundary moves closer to the main body of Fermi 
surface thus reducing the length of the cylinders. The number of electrons 
in these positions would be proportional to the circumference of the necks 
and hence be increased with the diameter of the area of contact. It is 
of interest to note that the zone distortion which causes the higher-than- 
ideal axial ratios has no effect on the trend in vol/atom as a function of 
e/a; whereas, as discussed earlier, the overlap of electrons across the 
{10-0} faces apparently causes an overall increase in the vol/atom (fig. 1). 

The axial ratio results in fig. 2 and the Aa/(a— Aa) plots in fig. 4 indicate 
that, in addition to increasing the number of electrons per atom, the 
solute elements also have an effect on the band structure of the alloys. 
It is possible to correlate these results with the parameter 7 defined by 
Cohen and Heine (1958) as a measure of the change in band structure 
caused by the addition of the solute element. The 7 parameters, in units 
of ev, for the B-sub-group elements which form ¢ phases with Ag are 
listed at the bottom of fig. 1. It is to be expected that similar trends in 
lattice spacings will be observed in those ¢ phases for which the y para- 
meters are nearly the same. Consider, for example, the alloys containing 
In, Sn and Sb from the second long Period and for which the size factors 
(fig. 1) are all favourable. The rate of change of the a’ spacings for the 
primary solid solutions in these systems are almost identical, whereas 
in the ¢ phases the rate of change of a (as shown by the plots of Aa/(a— Aa) 
in fig. 4) and of the decrease in c (fig. 1) and c/a (fig. 2) is much greater 
in Ag—In for which 7=5-7 than in Ag-Sn or Ag—Sb for which 7= 2-4 
and 2-9 respectively. If, on the other hand, alloys with solutes in the 
same B-sub-group are compared, it is found that Ag-Ga and Ag—h, 
which have similar 7 parameters, show similar rates of change of axial 
ratio and fractional distortion Aa/(a— Aa) as shown by the plots in figs. 2 
and 4 respectively. Again comparing solutes from the same B-sub-group. 
the changes in c/a and Aa/(a— Aa) with electron concentration are more 
rapid in Ag—Hg (y=5-9) than those in Ag—Cd (y=3-0). Ag—As (y= 3:7) 
forms the only exception to a consistent correlation between 7 parameters 
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and lattice spacing changes, since the changes in c/a and Aa/(a—Aa) in 
this system are more gradual than those in Ag-Sb (7 =2-9), even though 
the 7 parameter of the former is slightly greater than that of the latter. 
However, this discrepancy may not be too serious since the s—p excitation 
energies used to calculate 7 for the polyvalent solutes are known to be 
open to considerable error (Cohen and Heine 1958). 

Using the relationship for the energy of free electrons, EH=h?k?/2m, 
and that relating the electron concentration z with the radius & of the 
Fermi sphere, z=2Q. 47k? where Q is the atomic volume, an order of 
magnitude estimate of the band gap across a zone face can be made 
from a knowledge of the lattice parameters and the value of z at overlap. 
Assuming that the Fermi surface is very nearly spherical in Ag-based 
alloys and that z at overlap is given by the plots of Aa/(a—Aa) in fig. 4, 
sufficient information is available to evaluate the band gaps, AZo), 
across the {10-0} zone faces. A typical calculation for the Ag—Sn system 
for which z at overlap ~1-395 yields AH,,,.,~0-73ev. Further, since 
it has been shown that the band gaps divided by the energy of the 
electrons at the {10-0} and {00-2} zone faces are in the ratio 1:1-6 
(Goodenough 1953), the band gap across the {00-2} faces, AF io.) 
may also be obtained. Thus, for Ag-Sn AZ,..,,~1:32. Comparing 
this value with that for the band gap across the equivalent zone faces 
in the f.c.c. zone, for which Tibbs (1938) has calculated AH,,,,,~ 2-43 ev, 
it may once again be concluded that the addition of Sn (or any other 
B-sub-group solute) to Ag results in a decrease in the band gaps. 


ACKNOWLEDGMENTS 


We are grateful to Professor M. H. Cohen and to Dr. L. F. Vassamillet 
for many stimulating and informative discussions. We also wish to 
thank Mr. C. L. Ruffner for his valuable assistance with the experimental 
work. Following the submission of the paper, we have received several 
helpful suggestions from Professor Hume-Rothery, for which we should 
like to thank him. 


APPENDIX 
1. MATERIALS AND MrerHops 


All alloys were prepared from spectroscopically standardized materials 
(supplied by Messrs. Johnson, Matthey and Co., Ltd.). Weighed quan- 
tities of the pure metals were melted with frequent vigorous shaking 
under reduced pressure of helium in small, sealed, thin-walled tubes of 
fused quartz or “Vycor’; details of this technique have been published 
earlier (Raynor and Massalski 1955). Following previous experiences 
the synthesized compositions of the alloys were accepted without chAniteal 
analyses provided that no significant change in weight had occurred 
during casting. After casting all ingots were given a homogenizing 
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treatment (see table) in sealed pyrex tubes under a reduced atmosphere 
of helium, and later examined metallographically under polarized light 
for the presence of impurities and additional phases. X-ray powder 
Specimens of the alloys were requenched from the respective homogenizing 
temperatures. 

A 2 radian film camera with 0-5mm pin holes was used for the x-ray 
work, and the powder specimens of less than 325 mesh were held in 
03mm diameter capillaries. The film was used in the asymetrical 
position and measured with a vernier device (Klug 1940) which enabled 
line positions to be located to 0:02mm. The last ten high-angle lines 
in the diffraction patterns were used to calculate the lattice spacings by 
an extrapolation procedure involving successive approximations of the 
axial ratio (Massalski and King 1960b). All lattice spacings are reported 
in Angstrom units, based upon the following values assumed for CuKa 
radiation : 

CuKa, = 1:54050064 
CuKa, = 1:5443412 4. 


The temperature of the specimens during the exposure was 30°+1°c. 
The accuracy of the reported results is + 0-0002A4 for the a spacings, and 
+0-00034 for the c spacings. This accuracy of one part in 15000 is 
not absolute, but is indicative of the degree of reproducibility obtained 
in the authors’ laboratory. The significance of absolute and relative 
accuracy in the measurement of lattice parameters has been discussed 
recently by Parrish (1960). 


2. PRESENT RESULTS COMPARED WITH PREyIoUS Data 


A close comparison between the present results and the previous 
data available in the literature for silver-based ¢ phases is not possible, 
since most of the latter were obtained as a by-product of phase identi- 
fication studies and are often reported only to the third decimal place. 
The data reported by Pearson (1958) for several alloys in the systems 
Ag—Cd, Ag-In, Ag-Sn and Ag-—Sb, and by Hade and Hume-Rothery 
(1959) for Ag-As, are in general agreement with the present trends in 
lattice spacings, even though individual values may differ from those 
of the present work. Similar general agreement is found between the 
present results and Stenbeck’s data for Ag—Hg (Pearson 1958) even 
though only two values were reported. A single value reported by 
Rayson and Calvert (1958) for a Ag—Hg alloy is in complete disagreement 
both with the earlier work of Stenbeck and with the present results. A 
similar discrepancy between the results of Rayson and Calvert (1958) and 
those of other investigators has been reported previously for the lattice 
spacings of the ¢ phase in Au-Hg (Massalski and King 1960a). The 
values of both a and ¢ lattice spacings for a single Ag—Ga alloy reported 
by Hellner (Pearson 1958) are considerably greater than the present results 

for this system. Although it has been reported (Hume-Rothery and 
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Andrews 1959) from the results of high-temperature x-ray studies that 
the ¢ phase in Ag—Ga extends up to about 32at. %, in the present work 
it was found that this phase could not be retained, by quenching to room 
temperature, in alloys containing more than 26-5at.% Ga. 
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Free Precession Technique 
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THE idea that the inertial mass of a body may arise from a gravitational 
coupling with distant matter has led to suggestions that the concentration 
of matter near the centre of our galaxy may result in an anisotropy of 
inertia at the earth which could be detected experimentally. In particular, 
Cocconi and Salpeter (1958, 1960) have pointed out that such an anisotropy 
might cause shifts in energy levels of atoms and nuclei in magnetic fields. 
They suggest, for example, that in a nucleus which can be regarded as a 
single ps3). nucleon moving in an isotropic central potential the energies 
of the states with magnetic quantum numbers m,= + 3/2 would be increased 
slightly if the magnetic field were parallel to the direction of the centre 
of the Galaxy, and the energies of the states with m,;= +1/2 would be 
decreased by an equal amount; while if the magnetic field were perpen- 
dicular to the direction of the centre of the Galaxy the energies would be _ 
shifted in the opposite directions. From models of inertia like those 
discussed by Cocconi and Salpeter, as well as those referred to by Dicke 
(1959), one could expect that the ratio of the anisotropic part of the 
inertial mass of a body to the isotropic part might lie in the range from 
10-7 to 10-8; and in this case the suggested energy shifts would be large 
compared with the shifts which can be measured by nuclear resonance 
methods, being about 30Mc/sec even for the smaller anisotropy. In 
fact data from nuclear resonance experiments with crystals indicate 
that no energy shifts as large as this take place (Drever 1960). However, 
in view of the cosmological interest of the question, and the speculative 
nature of the calculations, it seems worth while to make the search for 
effects of this kind as sensitive as possible. Some relevant observations 
on lithium 7, a nucleus of spin /=3/2 with an unpaired pg. proton, are 
described here. 

Since this work was begun Hughes ef al. (1960) have reported a rather 
similar experiment done by nuclear magnetic resonance, which gave an 
upper limit to anisotropy of inertia of 10-°°; while Sherwin e¢ al. (1960) 
have obtained an upper limit of 5 x 10-16 from studies of resonance scatter- 
ing of gamma rays from *’Fe. The present experiment, which is more 
sensitive, has also given a negative result. 

The four energy levels of a nucleus of spin 3/2 in a magnetic field are 
normally taken to be equally spaced, and any modification of the spacing 
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by shifts like those suggested by Cocconi and Salpeter would split the 
single resonance obtained in a magnetic resonance experiment into a 
triplet. The minimum splitting detectable with simple resonance technique 
is usually limited by inhomogeneities in the magnetic field ; but in the 
present case the splitting, if small compared with the separation of the 
levels, is independent of the magnitude of the field and this restriction 
may be eased by working in a weak field. The magnetic field of the earth 
itself is very suitable, and this was the field used here. 


Fig. 1 
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Schematic diagram of apparatus. The switches, operated by relays, disconnect 
and short-circuit the amplifier while current is passed through the sample 
coil; and the Metrosil non-linear resistor limits the peak voltage generated 
during switch-off without affecting the nuclear signals. 


The experimental technique was based on that employed by Packard 
and Varian (1954) and Waters and Francis (1958) to observe free precession 
of protons in the earth’s field, and the arrangement is shown schematically 
in fig. 1. The lithium, in the form of a solution, is contained in a poly- 
thene bottle surrounded by a coil placed with its axis perpendicular to 
the direction of the earth’s field. Initially a direct current is passed 
through the coil to produce a field of about 200 gauss at right angles to 
the earth’s field. When this is switched off rapidly the resultant nuclear 
magnetic moment precesses about the earth’s field and an alternating 
e.m.f. is generated in the coil. Under suitable conditions with a single 
resonance, this signal decays exponentially with a time constant equal 
to the transverse relaxation time 7’, of the spin system. If, however, the 
resonance were split into a close triplet one would expect the signal to 
exhibit beats, corresponding to interference between oscillations at the 
three resonance frequencies which would be detected in a steady-state 
experiment. The results of a detailed analysis by Das and Saha (1955) 
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of the analogous situation of free precession in the presence of a weak 
electric quadrupole interaction indicate in fact that when the earth’s 
field is parallel to the direction of the centre of the Galaxy, and the splitting 
is at its maximum, the amplitude of the signal would vary as 


{2 +3 cos (27tAf)} exp (—#/T,), 


where Afis the frequency separation between each satellite and the central 
component of the triplet, and t is the time from the beginning of the free 
precession. As the earth rotates the splitting, and hence the modulation 
of the precession signal, would alter; and a maximum splitting too small 
to give rise to obvious beats might show as an apparent change in the 
decay rate with sidereal time. At the place where the experiment was 
performed, latitude 55°54'N, longitude 4°30’ W, the declination and 
inclination of the earth’s field are 71°N and 9° W respectively, and the 
direction of the field passes within 10° of the centre of the Galaxy at about 
5h sidereal time, when almost the maximum splitting would be expected. 

The frequency of precession of “Li in the earth’s field is only about 
800 cycles/sec, and in practice the signal was weak and it was necessary 
to use a large sample (2-5 1.) to obtain an adequate signal-to-noise ratio. 
The sample coil was situated in the open air in a position remote from 
steel-framed buildings and power cables; and to minimize pick-up of 
electromagnetic interference a second similar coil was mounted near the 
sample coil and connected in opposition with it. It was necessary to control 
the temperature of the lithium solution to maintain a constant relaxation 
time, and this was done by passing a suitable current through the coil 
in the intervals between measurements. A stirring device of flexible 
polythene was incorporated in the bottle, and the temperature was kept 
at (37+1)°c. Data on relaxation times (Bloembergen e¢ al. 1948) suggested 
that a solution of lithium nitrate would be a particularly suitable material 
to use as sample, and after preliminary tests with other compounds this 
was employed throughout. Nitrogen was bubbled through the solution 
before use to reduce the amount of dissolved oxygen, and the polythene 
bottle was hermetically sealed to prevent changes in oxygen content during 
the experiment. When care was taken to eliminate magnetic material 
from the vicinity of the sample (and some borosilicate glass proved appre- 
ciably ferromagnetic) the decay constant of the observed ‘Li signal was 
3:7 Sec. 

Two continuous 24-hour runs and numerous shorter ones were made 
(during May and July 1960), the precession signals being examined on an 
oscilloscope and recorded on film at intervals of 20 or 30min. A typical 
photograph is shown in fig. 2+. No beats of the type expected were seen, 
and no significant changes in the decay rate of the signal or its initial 
amplitude were detected. Some small changes in the shape of the envelope 
were in fact noticed but these were comparable with fluctuations due to 


+ Figure 2 is shown as a plate. 
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noise, were not reproduced from one run to another, and appeared to be 
correlated with slight changes in the temperature of the solution. The 
examination of the photographs was made by projecting them onto a 
drawing of the envelope expected for a pure exponential decay in the 
presence of noise. Comparision with a theoretical envelope for the case of 
a splitting for which Af = 0-04 cycles/sec (in which the first minimum occurs 
at t= 10sec) indicates that a slowly-varying splitting of this magnitude, 
which would arise from individual level shifts of 0-02 cycles/sec, would 
have been readily detectable. 

These results themselves do not completely rule out energy shifts of 
such a magnitude that the outer components of the triplet remained outside 
the pass-band of the amplifier for most of the sidereal day, and only came 
within it for short periods which happened to coincide with intervals 
between observations. A separate experiment was therefore performed in 
which the initial amplitudes of the precession signals due to the ‘Li 
nuclei and to the protons in the same lithium nitrate solution were com- 
pared. In doing this it was necessary to ensure that the decay of the 
magnetizing field was sufficiently rapid to give maximum signals from 
both types of nucleus, and to take into account the difference in the 
sensitivity of the apparatus at the two frequencies concerned (803 cycles/sec 
for 7Li and 2068 cycles/sec for protons). The experimental results agreed 
to within 5° with the ratio of the amplitudes calculated assuming detection 
of the whole “Li signal; and in disagreement with the ratio which would 
be expected if only the central component of a triplet were observed, which 
is smaller by a factor of 2/5. 

It may be remarked that if there were energy shifts large compared 
with the spacing of the levels in the earth’s field there would, in general, 
be no signal which remained at the expected “Li precession frequency 
throughout the sidereal day. 

It may be concluded that any shifts of the 7Li energy levels of the type 
suggested by Cocconi and Salpeter do not alter the spacing of the levels 
by more than 0-04 cycles/sec. If one applies the calculation of Cocconi 
and Salpeter directly to “Li, taking the mean kinetic energy of the pas). 
proton to be 10Mev, this result would correspond to an upper limit for 
the ratio of the anisotropic part of the inertial mass of a proton to the 
isotropic part of the order of 5x 10-8, 

It should be noted, however, that a negative result obtained in an 
experiment of this kind does not rule out the possibility of a much larger 
anisotropy in inertial mass which is masked by some other anisotropy. 
Epstein (1960) has recently suggested that an anisotropy in the potential 
energy of a nucleon might accompany an anisotropy in its mass and 
counteract it, and such a situation might well exist without being detected 
in an experiment like the present one. The experiment does, nevertheless, 
provide a check of the specific model proposed by Cocconi and Salpeter ; 
and in fact it is difficult to reconcile the shifts in energy levels suggested 


by this model with the very low upper limit for such shifts found experi- 
mentally. 
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ABSTRACT 


An account is given of palaeomagnetic measurements on borehole specimens 
from three dolerite sills of Carboniferous age. The magnetic directions were 
closely grouped for one sill but widely scattered for the other two, showing 
that contrary to the experience of previous workers, specimens from under- 
ground do not always yield satisfactory palacomagnetic results. 

The intensities of magnetization were very low at the margin of each sill. 
In one instance this phenomenon was studied in detail. It seems to be due to 
the fact that at the margins, the iron occurred mainly in the form of a non- 
magnetic carbonate, produced when the sill was intruded into limestones 

Some wider implications of the work are discussed. 


§ 1. INTRODUCTION 


THE present paper describes a series of magnetic measurements on material 
from deep borings through three dolerite sills in the Carboniferous system 
of Great Britain. The first boring was through a 400 foot sill at Rashiehill, 
Stirling, lying between depths 1353 feet and 1742 feet in beds of the Lime- 
stone Coal group; the second was through a 140 foot sill at Outmains, 
Ayrshire, between depths 1774 feet and 1915 feet in Upper Coal Measures ; 
and the third through a 40 foot sill at Plungar, near Nottingham, between 
3165 feet and 3205 feet, also in Upper Coal Measures. From the 40 foot 
sill complete fitting sections of 44in. diameter were available over most of 
the length. From the 140 foot sill, seven short sections 2in. diameter and 
from the 400 foot sill, 42 sections 44in. diameter were available. The 
depths of these pieces were known and the tops had been marked by the 
drillers. These markings will be referred to hereafter as the nominal tops. 
For the 400 foot-sill some of the overlying baked rocks were also available ; 
and in addition oriented samples were collected at Torpichen, near 
Bathgate, where the sill crops out. 

The original object of the work was to obtain an independent test of the 
palaeomagnetic directions obtained from other work on the Carboniferous 
system (Everitt and Belshé 1960). Measurements by Gough (1959) on 
the Pilansberg dykes of South Africa and by Jaeger and Joplin (1955) on 
Tasmanian sills had shown that the scatter of magnetic directions could 
be greatly reduced by taking material from underground rather than the 
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surface, presumably because it eliminated the effects of recent Weanhenas 
and lightning. The orientations of the borehole specimens were not known ; 
but for the present purpose this was not important, since the main point 
at doubt in the Carboniferous system was the magnetic inclination. 


Fig. 1 
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Magnetic directions in the three sills. 
O upward inclination. 
@ downward inclination. 
x declination relative to arbitrary arrow. 
(a) 400 foot sill, inclination; (b) shale baked by 400 foot sill, inclination; 


(c) 140 foot sill, inclination; (d) 40 foot sill, inclination; (e) 40 foot 
sill, declination relative to arbitrary arrow. 


Each sill was found to have very low magnetic intensities at the margins. 
Similar low intensities have been observed at the margins of the sills by 
Jaeger (private communication) and Hatherton (1954), and at the margins 
of dykes by Leng (1955) and myself (unpublished). Although Bruckshaw 
and Robertson (1949) and others have found examples of the opposite 

effect, that is an increase of intensity towards the margins of igneous 
bodies, the present phenomenon seems sufficiently widespread to be worth 
investigating. 


The palaeomagnetic directions are discussed in § 2 and the measurements 
on magnetic properties in § 3. 


§ 2. DrrEections or MAGNETIZATION 


The directions of magnetization of the original pieces of core, measured — 
on an astatic magnetometer are presented in fig. 1. For the 40 foot sill, 
there were complete fitting lengths extending over most of the core, so 
histograms are given for inclination and for azimuth relative to an 
arbitrary line on the core. For the other sills, with broken lengths, only 
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the inclinations could be measured. There is some uncertainty whether 
these inclinations are up or down, owing to possible mis-marking of the 
specimens, so they have been plotted without regard to sign, the measured 
upward and downward values being distinguished only by different symbols. 
Also included in fig. 1 are the results from a group of strongly magnetized 
baked pyritous shales from a depth 80 feet above the 400 foot sill; for these 
there was no possibility of mis-marking since the pieces fitted together. 
The numerical data from the histograms and measurements on surface 
specimens from the 400 foot sill are set out in Table 1. 


Table 1 
Standard 
. Type of | Numberof| Mean deviation | Deviation 
Sill Site é 5 aa : : 
material | specimens dip of single | of mean 
readings 
feet 
40 | Plungar_ | Borehole 14 24° down 12° 4° 
140 | Outmains | Borehole 7 42° down 16° 6° 
Borehole 42 16° down 5° 0-8° 
igneous 
- 4-7, | Borehole 12 0° 22, 6° 
400 | Rashiehill LOS) 
Surface 6 4° down 30° 12° 
igneous 


The only specimens that show really close grouping are those from the 
boring through the 400 foot sill; the results from the other two sills are 
widely scattered and must be discounted for palaeomagnetic purposes. 
Of the specimens from the 400 foot sill, the inclinations of 30 were nominally 
downards and 12 were nominally upwards. One conceivable explanation 
of this grouping is that some specimens had undergone self-reversal ; but 
in view of the well-known difficulty of marking broken lengths of cores 
reliably, I am inclined to think that the effect arises from mis-marking. 
The downward dip, being more frequent, is presumably the true one. The 
mean value is 16° down. “sed, 

The specimens from the outcrop of the 400 foot sill were very widely 
scattered, with a mean direction 216° east of the true north and 4° down. 
Taking 16° down as the best value for the magnetic inclination, and 216° 
east as the azimuth, these results fit in well with the mean value for the 
Carboniferous period as a whole of 200° east and 27° down (Everitt and 
Belshé 1960). 

We see therefore that for the 400 foot sill the scatter was greatly reduced 
by taking borehole specimens. For the other two sills, however, the 
scatter in the cores was large. Therefore sampling from underground 
does not always yield satisfactory palaeomagnetic results. On the 
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other hand, these results serve to emphasize that the existence of magnetic 
scatter is an interesting problem in itself, and that simple explanations 
of it in terms of recent weathering or lightning are inadequate. 


§ 3. InTENsITIkES OF MAGNETIZATION 


The intensities of natural magnetization for the three sills are plotted 
against depth in fig. 2(a), (b) and fig. 3(a). At the centres of the sills the 
intensities fell within the range 200-2000 x 10-*e.m.u./g, which is typical 
of dolerites; but at the margins, they were much lower. The greatest 
variation occurred in the 400 foot sill, where the least intensity was 
1:8 x 10-%e.m.u./g at the upper margin and the greatest was 5600 x 10~ at 
32 feet below the top. Most of the investigations were made on this sill ; 
the main results being set out in fig. 3 (a)-(f). 


Fig. 2 


InGIO® emus gm) 


Depth (feet) 


Ju (xIO° e.m.u./gm) 


Intensities of natural remanence. 
(a) 40 foot sill; (6) 140 foot sill. 


In principle there seem to be four general classes of explanation for the 
differences of intensity between the margins and centre of the sill: 


(1) A variation in the total iron content. 
(2) 


: A variation in the proportion of iron existing in a ferromagnetic 
orm. 


(3) A variation in the magnetic material giving rise to differences in its 
saturation intensity. 


(4) A variation in the degree of natu ramet 
7 ral remanent m t 
single magnetic material. agnetization of a 


These possibilities will be discussed in § 3.1. 
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3.1. Physical Cause of the Variations in Natural Remanence 


The total iron content for the 400 foot sill is plotted in fig. 3(b) was 
practically uniform throughout. These measurements, which are due to 
Dr. A. J. E. Welch, were obtained by X-ray fluorescence counts on flat 
discs 2-3 cm in diameter (Parrish 1956). Similar results were obtained for 
the other two sills. Thus possibility (1) above is eliminated. 


Fig. 3 
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Measurements on 400 foot: sill. 


(a) natural remanence; (6) total iron content; (c) saturation remanence; 
(d) ration of natural to saturation remanence; (e) saturation remanence 
after heating to 1000°o; (f) demagnetization by alternating fields. 


In fig. 3(c), the magnetic remanence (/,) after saturation in a field of 
3000 oersteds is plotted against depth for the 400 foot sill; while the ratio 
of the natural remanence (J) to this saturation remanence is plotted in 
fig. 3(d). EKach point on these diagrams represents the mean value for 
several small pieces cut from the original lengths of core. The ratio of 
Jy/Jz is approximately uniform, showing that the main variation in Jy 
does not arise from differences in the degree of magnetization. Thus 
possibility (4) above is eliminated. There are however one or two small 
variations in J,/J, which will be discussed in §3.3. Once again, similar 
results were obtained for the other two sills. 

To decide between the two remaining possible explanations of the low 
intensities, specimens from near the upper margin and the centre of the 
sill were powdered, and the magnetic fractions extracted. Since the 
material at the extreme margin was too fine grained to separate success- 
fully, the upper specimen was taken at a depth of 6 feet, where the natural 
intensity was 30 x 10-°, or about 3%, of the mean value across the centre 
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of the sill. The proportion of magnetic material proved to be obviously 
much lower at the margin than at the centre. In table 2 measurements of 
the intensity of saturation remanence for the native rock and the magnetic 
fraction are set out; together with estimates of the percentage magnetic 
fractions by weight, obtained from the ratios of the two intensities. 


¢ 


Table 2 


Remanence of Remanence of Proportion by weight 
Depth native rock magnetic fraction of magnetic material 
(e.m.u./g) | (e.m.u./g) (%) 


6 feet : . ipl 
95 feet . . 
220 feet 


The saturation remanence in the magnetic fraction from a depth of 6 feet, 
is about 8°, of the values at the centres of the sill. However, it was clear 
under a microscope that the fine-grained material from the margin was 
very heavily contaminated with quartz; so the estimate of its saturation 
remanence is merely a lower limit, and the estimate for the proportion of 
magnetic material is, similarly, an upper limit. 

The material has been examined petrologically by Dr. P. H. 8. Stubbs. 
Measurements on thin sections gave the proportion of opaque minerals 
as 5-10°% at the top margin and 10% in the centre of the sill. The main 
opaque mineral at the margin has been identified by Dr. A. P. Millman, 
using reflection microscopy, as ilmenite; which is, if pure, non-ferro- 
magnetic, 

Some of the ilmenite was separated by centrifuging the powdered rock 
in tetrabromoethane (density 2-96). A close examination of the finely 
powdered fraction showed once again that it contained a very small 
proportion of highly magnetic material. On heating, this gave exactly 
the same thermal demagnetization curve as the rock in bulk (see § 3.3). 
Thus the ilmenite did not contribute to the original magnetization. 

I conclude, therefore, that the intensity variation across the sill was 
mainly caused by differences in the proportion of iron existing in a magnetic 
form. Further measurements’ to investigate whether there were any 
differences in the magnetic material will be described in § 3.3. 


3.2. Geological Origin of the Variations in Intensity 


When the rock from the margins of the 400 foot sill was heated to a 
few hundred degrees centigrade, there was a very marked increase in 
intensity. Figure 3(e) shows that after treatment for half an hour at 
1000°c the increase was so great that the saturation remanence became 
practically uniform throughout the sill. | On the other hand, when the 
_ Magnetic fraction from the margins was heated separately, the intensity 
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was 1:4e.m.u./g before heating and 1-1 e.m.u./g afterwards. Therefore 
the increase must have been due to the breakdown of an initially 
non-magnetic compound. 

Two plausible possibilities are a hydroxide (or hydrated oxide such as 
chlorophaeite) and a carbonate. Microscopic examination by Dr. P. H.S. 
Stubbs disclosed the presence at the margins of an abundance of interstitial 
carbonate, which he considers to be derived from the overlying limestone, 
probably by replacement of the original pyroxene. 

A sample from the margin was powdered and divided into two parts, 
of which one was treated with dilute HCl to dissolve the carbonate. The 
two samples were heated at 900°c for a quarter of an hour in evacuated 
quartz tubes. On cooling, the untreated powder became highly magnetic, 
while the other was unaffected. 

A specimen from the margin was heated to 1000°c for séveral hours, 
powdered and separated with a small bar magnet. X-ray diffraction 
photographs of the magnetic fraction revealed the presence of magnetite 
and two materials, each with ilmenite structure but slightly different 
lattice spacings. The simplest explanation is that the newly produced 
magnetic material consists of intergrowths of magnetite and ilmenite, 
modified in lattice spacing by solid solution; and that this material 
co-exists with the original ilmenite phase. This observation rules out 
the only likely alternative explanation of the intensity increase, that is, 
that the original ilmenite was reduced to magnetite during the breakdown 
of the carbonate. 

I conelude that the proportion of magnetic material at the margins of 
the 400 foot sill is low because the iron mainly occurs in the form of a 
non-magnetic carbonate. This was probably produced when the sill 
was intruded into limestones. 


3.3. Properties of the Magnetic Material in the 400 foot Sill 


In fig. 4 heating curves are given for the saturation remanences at 
various depths in the 400 foot sill. Across most of the sill the curves 
were identical giving a single Curie temperature of 550°c. In the top 
6 feet there was a change in the shape of the curves, probably indicating 
the presence of a second magnetic material with a Curie temperature 
of 300°c. The magnetic material has not been studied in detail, since 
it was too fine-grained to separate easily. 

The opaque minerals were examined under the reflection microscope by 
Dr. P. H. S. Stubbs. As mentioned in §3.1 there was a high proportion 
of ilmenite at the margins, which was in the form of clusters of particles 
5. to 20u in diameter. In the rest of the sill there were large skeletal 
grains of magnetite, 0-3-0-8 mm across, containing ilmenite intergrowths, 
with some patches of secondary alteration to maghemite. The Curie 
temperature of 550°c is consistent with the presence of magnetite. 

The magnetic moments of the virgin rock and magnetic fractions were 
measured in fields up to 10000 oersteds on a vibration magnetometer 
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designed by Blackett and Sutton (1956). Figure 5 shows that at Se 
margins of the sill, there is. some material with ee eae : 
susceptibility (25 x 10-6 e.m.u./g in the rock as a whole) in 2 e e 
the ferromagnetic material. ‘This may well be the fee oe be 
at the margins, since ferrous carbonate (piderite) has a pease, . 
susceptibility of about 40-120 x 10-%e.m.u./g (Chapman an artels 
1940). 


Fig. 4 
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Heating of saturation remanences, 400 foot sill. 
x depth 0-5 feet. 
@ depth 50 feet. 
@ depth 95 feet. 
O depth 201 feet. 


There were systematic variations of the coercive properties across the 
sill. Figure 3(f) shows the effect of alternating field demagnetization 
on the saturation remanence, giving contours of the percentage reduction 
produced by peak fields of 100 and 200 oersteds. These curves were used 
rather than straightforward measurements of coercive force to eliminate 
the effect of the large background of paramagnetic material. They 
indicate a great increase of hardness (that is to say coercive force) 
the margins of the sill. This effect may be related to the decreas 
size at the margins, since it is well known (Néel 1955, p. 220) 
force increases with decreasing grain size. 

The ratio of natural to saturation remanence (fig. 3 (d)) has already been 
mentioned. There is a sharp increase at the extreme margins, which may, 
like the increase in coercive force, be related to the small grain size. At the 


towards 
e in grain 
that coercive 
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region of low natural intensity 90 feet below the top of the sill, the ratio 
of natural to saturation remanence was much lower. This probably arises 
from a secondary change in the natural magnetization: for there was no 
corresponding effect in the laboratory thermoremanence or saturation 
remanence, and heating experiments on three separate specimens revealed 
the existence of two components in the natural magnetization, each in the 
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Intensity versus field; specimens from margins of 400 foot sill. (a) curves for 
native rock; (b) curve for magnetic fraction. 


horizontal plane, differing in azimuth by about 120°. The origin of this 
effect is unknown. It occurs in the region of largest crystals, which was 
probably the last part of the sill to cool. uh 
These measurements therefore show that there were only minor variations 
in the magnetic properties of the material in the 400 foot sill; but major 
variations in the amount. Such variations as do occur in the magnetic 


2Y 
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properties are probably related to the grain size. No detailed study was 
made of the material from the other sills. 


§ 4. CONCLUSION 


The conclusions to this work may be summarized as follows: 

(1) The measurements on magnetic directions showed that contrary 
to previous experience, sampling from underground does not always yield 
satisfactory palaeomagnetic results. The only closely grouped directions 
were those from the boring through the 400 foot sill, which yielded a mean 
magnetic inclination of 16° down. Since the mean azimuth from the surface 
rocks was about 216° east, this result fits well with the figure of 200° east, 
27° down already obtained for the Carboniferous system (Everitt and 
Belshé 1960). . 

(2) The intensity of natural remanence was very low at the margins of 
each sill. A detailed study of the effect was made for the 400 foot sill. 
It appeared to be due to the fact that at the margins the iron was mainly 
in the form of a non-magnetic carbonate, produced when the sill was 
intruded into limestones. 

(3) The magnetic material was very uniform throughout the sill except 
within a few feet of the top. This effect is interesting primarily by contrast 
with other work, especially on lavas, in which major differences in Curie 
point have frequently been found between specimens only a few centimetres 
apart. Such differences may be due to weathering, or alternatively to the 
difference between the conditions of formation of sills and lavas. It would 
be of great interest to examine sections from borings through unweathered 
lavas. Any further measurements on the magnetic materials in sills 
would probably best be concentrated on ones which had suffered extensive 
differentiation. 
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The Relative Hardness of the Hard Directions in Diamond+ 


By Emren M. WriKs 
Clarendon Laboratory, Oxford 


[Received December 20, 1960} 


ABSTRACT 


A method is described of determining the relative hardness of the hard 
directions in diamond, using a micro-abrasion tester and diamond impregnated 
cutting wheels. Cube, dodecahedron and octahedron planes have been 
studied and an order of hardness established between the various directions 
on these planes. This order of hardness should also apply for ordinary 
diamond polishing when loose diamond powder is used on a flat scaife, 
although in this case the differences in hardness are very much greater. 


§ 1. INTRODUCTION 


THE anisotropy of the hardness of diamond is very relevant to the shaping 
of diamond tools and to their use for working other hard materials. 
Diamond polishers have known for a long time that the dodecahedron and 
cube planes can only be polished in certain directions, and that the octa- 
hedron plane is extremely resistant to abrasion in all directions. Quanti- 
tative measurements of the resistance of diamond to abrasion by a cast- 
iron wheel charged with diamond powder, as is used in conventional 
polishing techniques (see for example, Grodzinski, 1953) have been made 
by Tolkowsky (1920), Kraus and Slawson (1939), Denning (1953, 1955, 
1957), Wilks and Wilks (1954, 1959), and Hukao (1955). As a result a 
considerable amount of information is now available regarding the vari- 
ations in the hardness in different directions and on different planes. 
Although there is general agreement as to which directions are most easily 
polished, there has been some confusion as to which are the hard directions, 
and it is this problem that we now discuss. 

The cube plane of diamond is very resistant to abrasion in the four 
directions at 45° to the crystal axes, the dodecahedron plane in the two 
directions at 90° to the axis, and the octahedron plane in all directions. 
On the latter plane, all the directions are not equally hard, and there has 
been some discussion as to which are the hardest (Denning 1953). Using 
conventional polishing techniques it is very difficult to abrade diamond at 
all in these hard directions, and there has hitherto been little information 
as to their relative hardness. Although there seems little prospect of 
making such measurements using the technique of a cast-iron scaife and 
loose diamond powder, it is now possible to compare the hard and soft 
directions on a cube face of diamond by using diamond impregnated 
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grinding wheels (Wilks and Wilks 1959). The variations in hardness so 
observed are much less than those observed using conventional loose 
powder techniques, but as we discuss later the same order of hardness 
should be observed in both measurements. We have therefore applied 
this technique to compare the hardness in the hard directions of the 
three principal planes : cube, dodecahedron and octahedron. 


§ 2, ExPpERIMENTAL METHOD 


As with the earlier measurements (Wilks and Wilks 1959) a Grodzinski 
micro-abrasion tester was used with specially prepared diamond impreg- 
nated wheels of approximately lin. diameter mounted on accurately 
centred spindles. The wheels rotated at a speed of 10000r.p.m. and were 
brought into contact with the diamond under test under a load of 200g 
for 20secs; olive oil being used as a lubricant. Smooth cuts were pro- 
duced whose depths were subsequently determined by multiple beam 
interferometry ; the inverse of the depth being taken as a measure of the 
abrasion hardness. 


Two octahedron diamonds of fair quality and of size about 2 carats were 
chosen for the tests. On each stone we were able to select octahedron 
planes with very few surface markings. Cube and dodecahedron faces 
were cut and polished on the stones, and checked that they were aligned to 
within 30min of the lattice plane (when it was necessary to repolish these 
planes between two series of tests, this alignment was maintained). 

In order to compare the depths of the cuts made in the hard directions 
on these three planes, it is necessary to allow for any wear of the wheel. 
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Cuts were therefore also made in the soft direction on the dodecahedron 
plane. The figure shows the layout of a typical set of 14 consecutive 
abrasions, cuts 1, 2, 5, 6, 9, 10, 13 and 14 serving as calibration cuts. 
Finally for completeness cuts were also made in the principal directions of 
the three main planes. 


$3. RESULTS 


Table 1 shows the average of the depths of each pair of cuts made in 
the various directions in units of \/2. The two cuts of a pair had equal 
depths to within 10%, which as before (Wilks and Wilks 1954) sets the 
limit of accuracy of the experiment. Each line in the table represents a 
series of cuts made with the same impregnated wheel on the same diamond. 


Table 1. Average depth of cut in units of X/2 


Plane Dodecahedron Cube Octahedron 
(011) (001) (114) 

Direction x tig 7 ss 

of abrasion [100] [O11] [100] [110] [211] [211] 


Diamond Wheel 


If the depth of cut on the dodecahedron plane in the soft direction is now 
taken as 10-0 \/2, the depths of the other cuts can be normalized to this 
value, and the average of these results is shown in table 2. 


Table 2 


Plane and direction 


Plane Dodecahedron Cube Octahedron 
(011) (001) _ Pg ah ae a 

Direction [100] [O11] [100] [110] [211] [2TT} 

Depth of cut A/2 | 10-0 1:8 8:5 1-1 3°9 2-1 
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The resistance to abrasion is the inverse of the depth of cut made by the 
abrading wheel. If therefore we take the resistance to abrasion of the 
soft direction on the dodedcahedron plane as 1-00 the results in table .2 


can be presented in order of increasing hardness, as shown in table 3. 


Table 3. The resistance to abrasion of various planes and directions 
presented in order of increasing hardness 


Plane (O11) (001) (111) (111) (011) (001) 
Direction [100] [100] [211] [211] [O11] [110] 
Hardness value 1:0 1-2 2-5! 4:5 5:5 8:5 


§ 4. Discussion 


The present experiments indicate the relative hardness between the 
various directions in diamond as measured with diamond impregnated 
wheels; now we wish to relate these results to the resistance of diamond 
to abrasion determined by conventional polishing techniques. It has 
previously been shown that comparable results are obtained from experi- 
ments made either with a flat-scaife charged with loose diamond powder 
or with the cast-iron wheel of a micro-abrasion tester also charged with 
loose diamond powder, e.g. the variation of the rate of abrasion with speed 
and load, and the effects of orientation (Wilks and Wilks 1954, 1959). In 
the present experiments the abrading diamond dust was bonded into the 
wheel, not as in normal polishing where the wheel is charged with loose 
powder in olive oil. Previous experiments have indicated that the use of 
loose powder magnifies the basic hardness differences which in themselves 
are quite small (Wilks and Wilks 1959). This view is confirmed by the 
results in table 1 which show that the relative hardness between the soft 
directions in the cube and dodecahedron planes is in agreement with what 
we have previously observed using loose powder. Thus although the 
hardness differences when loose powder are used are very much greater, 
the order of hardness for the various directions is the same in both cases. 
Thus the order of hardness shown in table 3 for impregnated wheels should 
also apply for ordinary diamond polishing when loose powder is used on a 
flat grinding scaife. In particular the hardest directions are those on the 
cube plane at 45° to the axes. 

As we mentioned above, there has been some discussion as to which are 
the relative hard and soft directions on the octahedron plane. Tolkowsky 
was unable to grind this plane at all, but by working on planes close to it, 
he deduced that the hardest direction was that away from the dodecahedron 
plane. Similar results were later obtained by Slawson and Kohn (1950), 
but Denning (1953) by actually grinding a plane parallel to the octahedron 
paar me Beau, oe the easiest, and the opposite direction the 
Seige eens re 8 confirm that the hardest direction is away from the 
ears ae Bs ethene with Tolkowsky, and with Slawson and 
ee ane ae ed that k oth directions on the octahedron plane are 

ard directions on the cube and dodecahedron planes. 
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Drrect electron microscope examination of thin foils prepared from 
neutron bombarded copper has revealed the presence of large numbers of 
small defects and dislocation loops (Silcox and Hirsch 1959). Similar 
effects have been observed in copper bombarded with «-particles, and in 
the latter case the occluded helium was found to modify the annealing 
behaviour (Barnes and Mazey 1960). While current theories of radiation 
damage require that primary knock-ons should be responsible for most of 
the damage remaining in a metal after bombardment, much interest has 
also been shown in focusing collisions which it is thought can transmit 
interstitial damage over quite large distances at energies less than 1 kev. 
Neutron bombardment experiments are difficult to analyse because of 
uncertainties in the neutron energy spectrum and specimen temperature, 
while the use of high energy ions or «-particles does not permit one to 
distinguish between damage occurring at high energies and that occurring 
after energy degradation. Itseems desirable, therefore, to have information 
on the damage produced by particles of much lower energies. However, 
most of the published work on low-energy ion bombardment is concerned 
with sputtering yield measurements and little attention has been paid to 
sub-surface damage produced in the crystal lattice. This note is a pre- 
liminary report on the nature of this sub-surface damage and its relation 
to radiation damage produced by more conventional techniques. 

Thin films of annealed gold 99-998°% pure, prepared by electro-polishing 
in a chromic acid—acetic acid bath, were mounted for electron microscope 
examination. The mounted foils were then bombarded in a specially 
designed ion gun. | 

The argon ions (A*) used in these experiments had an energy distribu- 
tion with a sharp high energy cut-off at 75ev and a half peak width of 
5ev. (The maximum energy a 75ev argon ion can transfer to a gold atom 
is 42ev in a held-on collision.) The angular spread of the beam was 
estimated to be less than 8°. The measured total current flowing to the 
specimen is due to the ions striking the surface and also to the secondary 
electrons leaving it. Assuming a secondary electron coefficient of 1 electron 
per incident ion, the total dose delivered to each specimen was 3-5 x 10" 
ionscm-2 in two hours. The background pressure of argon was 4 x 10~° 
mm Hg, and of contamination 1 x 10-*mmHg. The temperature of the 
specimen could be controlled during the experiment. 
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Figure 1} shows an area that was shielded from the beam on a specimen 
bombarded at +22°c. The structure is typical of a fully annealed foil. 
There are no dislocation loops or other defects, but there is a fine polishing 
structure on the surface. 

Figure 2 shows a bombarded area from the same foil. The contrast 
may be interpreted as being due to small planar defects, probably disloca- 
tion loops, which show up as short straight lines, sharp if the defect is 
looked at edge on. These are surrounded by their strain fields, a pair of 
dark or light crescents one on either side of the line (fig. 3). Similar contrast 
has been observed around quenched in dislocation loops in aluminium 
(Silcox 1960, private communication), and also around G.P.II platelets 
in an aluminium copper age-hardening alloy (Nicholson and Nutting 1958), 
although in neither case was the strain field as marked. 


Proposed interpretation of the contrast in fig. 2. 


The crystallographic orientation of the foil was found from a same area 
electron diffraction photograph, and where the defects line up head to tail 
to form a string, they do so along (111) directions. There is only one such 
direction near the plane of the foil in fig. 2. As well as the larger defects of 
100A diameter there are also smaller defects randomly distributed going 
down to the limit of resolution, less than 204. 

Figure 4 shows another grain with a different orientation in the same 
foil. Here the general defect size is smaller, and the distribution more 
random. The detailed effects produced by a particular bombardment were 
found to be very dependent on orientation. 

Figure 5 shows a grain from a second gold foil bombarded to the same 
dose at a lower temperature (—32°c) and subsequently examined at room 
temperature. The concentration of defects is much higher and many of 
the larger ones can be identified as dislocation loops. Also present in this 
grain are irregular strings of loops and short heavily jogged lengths of 


dislocation. The background of very small defects is more pronounced 
than in the first specimen. 


+ Figures 1, 2, 4 and 5 are shown as plates. 
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In interpreting the second specimen it is well to remember that it does 
not represent the structure as it existed at the temperature of bombard- 
ment, but only the structure after it has subsequently been allowed to 
anneal at room temperature. Consequently defects frozen in at the lower 
temperature will be able to move subsequently, and this may well account 
for the dislocations present being heavily jogged. 

There is a marked decrease in defect density towards the thinner outside 
edge of all foils examined, and this can also be seen in fig. 5. The thickest 
gold film in which detail is visible is about 15004; the foils shown here are 
all less than this and tapering. Although there were marked differences 
in the sizes and density of defects in different grains, in no specimen was 
any grain boundary effect detectable, either preferred growth or denudation 
in the boundary region. 

The efficiency of point-defect production was estimated by comparing 
a figure for the integrated planar defect area per unit area of specimen with 
an estimated value for the incident ion current. For a grain in the foil 
bombarded at +22°c, the efficiency was of the order of 5 x 10-4 point 
defects retained per incident ion. 

The penetration of 75 ev argon ions into gold can be estimated from the 
effective collision radius as used by Silsbee (1957). He derives an expression 
for the effective atom radius in an inter-atomic collision as a function of 
the energy transferred to the struck atom. For gold for 40 ev transferred 
the effective atom radius is about 0:55r, for 5ev transferred 0:7 r, where 2r 
is the equilibrium distances of closest approach in the gold lattice. These 
contractions are not sufficient to allow the argon ions to penetrate beyond 
the second layer of atoms. 

It is possible that there might be some diffusion-assisted penetration 
of argon. The observed efficiency of point defect production is low and 
does not necessarily rule out this possibility. However, Tucker and Norton 
(1960) observed appreciable diffusion rates only above 100°c whereas 
the effects observed here are in the —30°c to + 20°c range where diffusion 
rates will be very muchslower. Inno case has anything been seen resembling 
the small bubbles of helium gas reported in copper after «-particle irradia- 
tion (Barnes and Mazey 1960). 

If not due to argon atoms, the damage must be due to interstitials or 
vacancies, or possibly both. The temperatures considered lie above the 
recovery range for interstitial atom damage as determined from resis- 
tivity measurements and at the bottom end of the temperature range for 
vacancy migration (e.g. Seeger 1958, Seitz 1959). Interstitial atoms would 
be expected to diffuse very rapidly at all temperatures in the range 
investigated, while vacancies could diffuse only slowly at the lower 
temperature but considerably faster at the higher. 

The marked difference in the scale of the damage seen at these two 
temperatures suggests that vacancies are primarily responsible. However, 
vacancy migration could be merely serving to modify a structure produced 
initially by interstitial atoms. 
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If interstitial atoms are the cause of the damage, they may have either 
diffused in from the surface or they may have been injected to a depth of 
many atomic planes by a focused collision sequence of the type proposed 
originally by Silsbee (Silsbee 1957, Nelson and Thompson 1961, Vineyard 
et al. 1959 and 1960). Interstitials created by this second process are less 
likely to be lost to the surface before being captured by a sink inside the 
material, and therefore will be more effective in producing sub-surface 
damage. Further work will be needed to clarify the situation. 

Whatever the detailed processes involved, the results show that even 
in this low-energy range damage may be produced in gold in appreciable 
quantity by ion bombardment. 
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REVIEWS OF BOOKS 


Radio Astronomy. By F. Granam Smitu. (Penguin Books, 1960.) [Pp. 265.] 

Price 7s. 6d. 

Dr. Graham Smith’s book is an easily readable account of radio astronomy 
which can be understood without any specialized radio or technical knowledge. 
Details of radio receivers and radio telescopes are in fact left to closing chapters. 
The book covers the whole field of this rapidly advancing subject. Although 
Dr. Smith initially defines radio astronomy as “‘ the study of heavenly bodies 
by the reception of the radio waves they emit ”’ he does include descriptions of 
radio echo investigations or radar astronomy. A more inclusive definition of 
radio astronomy as “ the study of heavenly bodies by radio methods ” would 
seem preferable. Dr. Smith also discusses briefly the contribution to informa- 
tion about the ionosphere which has been derived through radio astronomy 
and by observations of earth-satellites. A final chapter outlines the main 
radio-astronomical observatories throughout the world. The book is well 
illustrated with diagrams and Plates. 

Dr. Smith enlightens the discussion with interesting commentary which 
helps one to view the development of the subject in proper perspective. Anyone 
with general scientific interests will appreciate the availability at such a reason- 
able price of a comprehensive, non-technical, and up-to-date picture of radio 
astronomy. : qe tepet a i 


Advances in Electronics and Electron Physics, Vol. XII: Photo-Electronic 
Image Devices. Editor: L. Marron. Volume editors: J. D. McGee and 
W. L. Wilcock. (Aeademic Press, 1960.) [Pp. 397.] Price $12.00. 

THis volume contains the papers contributed at a Symposium on Image 

Tubes and Related Devices organized by Professor J. D. McGee and held at 

Imperial College, London University, in September 1958. This publication is 

a most welcome inclusion of conference proceedings as part of a continuing 

series of volumes that usually are purchased as a series by scientific libraries. 

In this way the records of meetings on specialized topics are available for 

reference later when the field suddenly expands due to universal application 

of the original subject matter. This may happen in the present case when 
image conversion tubes have been improved beyond the technical limits 
discussed in these proceedings and are manufactured in quantity. The 
principles involved may be described as applied electron physics and the 
content is connected directly with the general theme of this series. Astronomical 
applications are discussed in the first three and two subsequent papers, out of 
the 30 contributions and x-ray uses in the last four. In between, various 
physical topics are presented in relation to image tubes such as signal-to-noise 
ratio, secondary electron emission, field emission and contrast perception. 

Intensification of the tracks of particles traversing a scintillation chamber is 

discussed as a possible application in nuclear physics. The only major omission 

is the area of solid state image intensifiers where a photocathode and phosphor 

are in contact and the applied field is across these two layers alone without a 

vacuum path for the emitted photo-electrons. No Russian papers are included 

despite extensive work in this general field reported previously. Hach paper is 
followed by the discussion from the floor which in many cases clarifies ambi- 

guities or omissions in the contribution. The experiment of publishing a 

symposium proceedings in this series has been successful and we hope will be 

repeated with other appropriate topics in electronics and electron Pe 
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Stochastic Processes: Problems and Solutions. By Lasos TakAcs. (Methuen 

& Co. Ltd., 1960.) [Pp. 137.] Price 18s. net. 

Tus unusual little volume is built around a collection of ninety-four problems, 
but it is rather more than a problem book. Each of the first three chapters 
(Markov chains, Markov processes, Non-Markovian processes) begins with a 
rapid but comprehensive sketch of basic theory, the problems then follow in 
well-graded sequence, and solutions are outlined in a fourth chapter. 

As a result of this arrangement, the book is self-contained to the extent 
that “all probems can be solved by using the listed theorems only”. How- 
ever, one would certainly be well advised to read a certain amount in stochastic 
processes beforehand, or to be prepared to do so concurrently, using Takacs’ 
book as a guide. The problems themselves are of a very reasonable standard: 
all worth while, some of them challenging, but, thanks in part to the systematic 
ordering of the examples, only a few seem really difficult. 

Topics covered in the theory include some spectral theory of transition 
matrices, the classification of Markov chains; diffusion, stationary, recurrent 
and renewal processes. The initial problems of each chapter are fairly direct 
illustrations of theory; then follow applications to queueing, telephone traffic, 
cascade processes, counter theory, statistical prediction, and some of the simpler 
statistical models of physics and astronomy. 

A little effort spent on Professor Takacs’ book will yield virtually nothing, 
but the result of a moderate effort should be quite unusually satisfying and 
rewarding. P. W. 


Paramagnetic Resonance in Solids. By Witu1am Low. Solid State Physics, 
Supplement 2. (New York and London : Academic Press, 1960.) [Pp. 212.] 
Price $7.50. 


PARAMAGNETIC resonance spectroscopy has become an important branch of 
Solid State Physics. The present book concentrates on a particular part of 
this field which has been intensively investigated during the past ten years, 
namely, paramagnetic resonance from transition group ions in magnetically 
dilute non-conducting crystals. The author has made many contributions in 
this area and he treats the subject clearly and extensively. After a short 
introduction, there is a useful chapter (67 pages) on the theory of crystal fields 
in these substances, in which are collected together tables of matrix elements, ~ 
group representations and operator equivalents together with some examples 
of how they are applied. This is followed by a chapter on the experimental 
results for these crystals (73 pages) which are extensively described with many 
references. There is some discussion of covalent bonding and also of optical 
absorption spectra in these compounds. The remaining three chapters (45 
pages) deal, rather briefly, with some of the many other aspects of paramagnetic 
- resonance in non-conducting solids, including colour centres, relaxation times 

masers and experimental techniques. Although these topics are given short 
treatment, and some others are omitted altogether, e.g. conductors, the book 
can be recommended as a good survey of work on transistion group ions, 


J. O. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents.] 


P. DELAVIGNETTE and S. AMELINCKX Phil. Mag. Ser. 8, Vol. 6, Pl. 83. 
Fig. 1 


(¢) 


Hexagonal prismatic loops in antimony telluride. (a) Very rapidly cooled 
(air quenched); the loops are rounded; () moderately slowly cooled: 
the loops are hexagonal; (¢) slowly cooled (~600°c/4 h) the loops form 
elongated hexagons. Notice the interaction between loops and glissile 


dislocations in the c-plane at the points marked A. 


P. DELAVIGNETTE and S. AMELINCKX Phil. Mag. Ser. 8, Vol. 6, Pl. 84. 
Fig. 2 


(0) 
Concentric pairs of elongated large loops in very slowly cooled specimens 
(~600°C/24 h). The two loops systematically differ 60° in orientations. 
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Electron micrograph of a Ni-20°% Co single crystal deformed at 300°K in stage I. 
Shear-strain interval Ae: 0-0-0612. 


Fig. 3 


Electron micrograph of a Ni-20% Co single crystal deformed at 90°K in stage I. 
Shear-strain interval Ae: 0-0-0724. 
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Fully annealed gold containing a grain Gold foil bombarded at 22°c. Disloca- 
boundary. The diffuse black lines tion loops lined up along the [111] 
are Bragg extinction contours. direction. . 


Fig. 5 


aha Bee in the same foil as figs. 1 Gold foil bombarded at —32°c. Disl 
aay i Bee a different arrange- tion loops and jogged dislocate 
. well as fine scale damage. 


